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Low- dimensional systems are an important field of current theoreti- 
cal and experimental research. Theoretically, the role of dimension- 
ality has been recognized for many years and dramatic predictions 
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have been made that still await experimental confirmation or are 
currently under study. Recent technological developments provide 
many possible realizations of effectively one-dimensional systems. 
These devices promise to give us access to a new range of phenom- 
ena. It is therefore very interesting to develop theoretical methods 
specific for such systems to model their behavior and calculate the 
correlators of the resulting theory. Incidentally, one such method 
exists and is known as Bosonization. It can be applied to one- 
dimensional systems and effectively describes low energy excita- 
tions in a universal way. It was developed in the 1970's when 
one-dimensional physics was viewed as a toy model for higher di- 
mensional physics. We use the example of a correlator known as the 
Emptiness Formation Probability to show that Bosonization fails 
to describe some long range correlators corresponding to large dis- 
turbances (the EFP measures the probability for the ground state 
of the system to develop a region without particles). We trace this 
failure to the fact that Bosonization is constructed as a linear ap- 
proximation of the full theory and we set up to develop a collective 
description with the required non-linearity. The resulting scheme 
is essentially a Hydrodynamic paradigm for quantum systems. We 
show how to construct such a hydrodynamic description for a vari- 
ety of exactly integrable models and illustrate how it can be used 
to make new predictions. For the special case of the spin-1/2 XY 
model we take advantage of the structure of the model to express 
the EFP as a determinant of a very special type of matrix, known 
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as Toeplitz Matrix. We use the theory of Toeplitz determinants 
to calculate the asymptotic behavior of the EFP in the XY model 
and discuss its relation with the criticality of the theory. Finally, 
we analyze the behavior of a charged particle in a two-dimensional 
medium filled with point-like magnetic vortices. 
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Preface 



This thesis is based on the work I performed during my graduate studies. 
Chapters [2] and [3] are devoted to the calculation of a correlator known as 
Emptiness Formation Probability for the anisotropic XY model. This analysis 
is based on the work published in [f] and [f]. In this chapter we also consider 
the effect on the block entanglement of the factorized ground state wavefunc- 
tion on a line of the phase-diagram of the model [§] . Chapter [6] describes how to 
develop a two-fluid description for a system of spin-1/2 electrons in one dimen- 
sion. This is still a work in progress [*] the ultimate goal of which is to address 
the problem of calculating the corrections to the exact spin-charge separation 
of the Luttinger Liquid model. Chapter [7] describes a two-dimensional prob- 
lem I studied with the help and the advice of Prof.A.S. Goldhaber [o]. We 
consider the Aharonov-Bohm effect for a scalar electron entering a medium 
filled with point-like magnetic vortices pinned to the sites of a square lattice 
and we consider the effect of such a configuration in the electron wavefunction. 
In Appendix [D] we explicitly construct the conserved currents of gradient-less 
hydrodynamic theories. While the existence of two infinite series of conserved 
quantities is known for these systems, we have not found such an explicit con- 
struction in the literature. This work is still unpublished, as the fundamental 
origin of these conserved quantities remains unclear (we address this issue in 
Appendix [D] but our results are so far inconclusive). 
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Chapter 1 
Introduction 



Dimensionality has an extremely important role in determining the physics 
of a system. Since the 1960's, one-dimensional physics stopped being just 
a theoretical toy, when excitations with unusual quantum numbers were ob- 
served in polyacetylene structures p]. Later on, in the 1980's, the discovery 
of novel effects in two-dimensional systems such as the Quantum Hall Effect, 
high-temperature superconductivity, and others, brought increasing attention 
to low-dimensional physics. 

In recent years, zero- dimensional ( "Quantum Dots") and one- dimensional 
( "Quantum Wires") systems have been implemented in laboratories as well, 
by effectively confining the electrons in fewer dimensions, with the motion in 
the transverse directions quantized so that at sufficiently low temperatures 
they are made energetically prohibited. 

From a theoretical point of view, one- dimensional models are particularly 
interesting, because, in some sense, all 1-D theories are strongly interact- 
ing, since particles cannot avoid each other in their motion (all collisions are 
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'head-on'). Moreover, in one dimension the distinction between a real scatter- 
ing event and the effect of quantum statistics is somewhat arbitrary, since it 
is impossible to exchange two particles without having them interacting. 

In addition to these remarkable aspects, there is something special about 
one-dimensional theories. The special symmetry of having just one spatial 
dimension together with the temporal one renders one-dimensional models 
particularly appealing and 'easy' to address, while preserving, and sometimes 
generating new, interesting non-trivial situations. For example, a number of 
theories are known to be exactly solvable in one dimension, quantum field 
theory methods are especially powerful and direct, and even the study of 
gravity is so simple in 1 + 1 dimensions that many attempts toward quantum 
gravity are performed in this lower dimensionality. 

In the 1970's, progress in quantum field theory was always parallel to the 
study of the corresponding 1-D theory, since the latter was often easier to 
tackle and provided important clues on the general structure of the methods. 
These efforts resulted in the development of the field theory method that we 
call "bosonization" and in the concept of the "Luttinger Liquid". 

1.1 Bosonization, linearization of the spectrum, 
and Luttinger Liquid concept 

One of the most successful methods to study one- dimensional systems is known 
as "bosonization" [15] which effectively describes low energy excitations in 1-D. 
One of its biggest advantages is its universality, i.e. the fact that the theory 
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has only two free parameters to be determined microscopically, when possible, 
and then the structure of the model is always the same for every system. 

This universality originates in the fact that the model assumes a linear 
spectrum, hence its applicability only at energy close to the Fermi points, 
and this assumption results in a quadratic Lagrangian, even including the 
interactions. 

Physically, the bosonization method provides a very powerful description 
of the concept of "Luttinger Liquid" (LL) [16] that is commonly employed to 
describe electrons in one dimension. This is the one-dimensional equivalent of 
the Fermi Liquid concept, but the fact that in one dimension the Fermi surface 
collapses to just two points has profound consequences on the physics of the 
system. Bosonization captures this physics accurately and allows for an easy 
access to the calculation of many correlation functions, often on very general 
grounds. 

The main limitation of the Luttinger Liquid description and of bosoniza- 
tion is the aforementioned assumption of a linear spectrum, which guarantees 
its universality. This approximation is extremely reasonable when one is inter- 
ested in low temperature physics, when only low energy excitations, close to 
the Fermi points, can contribute. But we would like to argue that there exist 
some important problems in one-dimensional physics where the non-linearity 
of the spectrum is essential and non- avoidable. In these cases one needs to go 
beyond the LL model. 

For instance, one of the problems we are going to address is the afore- 
mentioned prediction of spin-charge separation. It is a standard result of the 
Luttinger Liquid description that electrons in a one-dimensional system exist 
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in the form of spin and charge density waves whose dynamics is decoupled at 
low energies. As the LL model relies on the assumption of a linear spectrum, 
one can consider the effect of non-linearity of the spectrum on this prediction 
and this would couple spin and charge degrees of freedom. Any curvature 
in the spectrum would mix spin and charge degrees of freedom and if one 
looks at sufficiently high energy excitations the effect of this coupling would 
be measurable. Unfortunately, a perturbative treatment of the curvature of 
the spectrum in the bosonized theory generates divergences in the calculation 
of physical quantities and is notoriously difficult to perform. 

The bosonization description of a system is a collective description in which 
density fluctuations are the effective degrees of freedom considered. To over- 
come the assumption of linearity, we develop a collective description which 
retains all the characteristics of the spectrum, while still relying on just a few 
macroscopical fields. 

1.2 The hydro dynamic approach, integrable 
models, and Bethe Ansatz 

In fact, the collective description of quantum systems was suggested long ago 
by Landau [T7j. His method is essentially a hydrodynamic description in which 
the system is described just by its local density p and velocity v. These 
two fields obey a continuity equation and an Euler equation that depends on 
the original spectrum of the theory and encodes it into the dynamics of the 
hydrodynamic parameters. 



i 



We are going to apply this hydro dynamic approach to one-dimensional 
systems at zero-temperatures and show how to derive the hydrodynamic de- 
scription of different one- dimensional systems |19j . 

We will pay special attention to integrable systems. These theories are 
exactly solvable and have an infinite series of conserved quantities. Integrable 
models are very often good approximations for physical systems and an im- 
portant test field for theoretical methods, since their high degree of symmetry 
allows for very controlled calculations. 

Results for integrable models are usually derived within the framework of 
the Bethe Ansatz technique [201 [2T1 [T8] . The Bethe Ansatz provides a way to 
construct the wavefunctions of the system and to derive many thermodynamic 
quantities. The knowledge of the wavefunction, though, is implicit and this 
makes it very hard to calculate the correlation functions of the model. Exact 
expressions for the correlators are available in terms of determinants of oper- 
ator valued matrices [18], but these expressions are quite convoluted and are 
often not of practical use for the evaluation of physical quantities. 

We are going to show how the hydrodynamic approach derives directly 
from the Bethe Ansatz formalism. Therefore, physically relevant correlators 
involving the density and velocity of the fluid can be more easily evaluated 
within our approach. 

The ability of bosonization to calculate correlators for integrable models 
has been one of the great advantage of this method in conjunction with the 
Bethe Ansatz solution which provides the parameters of the theory. Clearly, 
these correlators are correct up to the linear spectrum approximation. The 
hydrodynamic approach overcomes this limitation at the cost of dealing with 
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a more complicated collective description, since the hydrodynamic equations 
are intrinsically non-linear. 

1.3 The Emptiness Formation Probability 

Throughout most of this work we will concentrate on a particular correlation 
function known as "Emptiness Formation Probability" (EFP). This correlator 
measures the probability that a region of the system is depleted of particles 
[TS] . It is clear that a bosonization approach is not suitable for the calculation 
of the EFP, since particles from everywhere in the spectrum are needed to 
empty a region of space and the curvature of the spectrum cannot be neglected. 

The EFP will be a useful example to show the advantage of the hydro- 
dynamic approach over bosonization when contributions far from the Fermi 
points are involved. But the importance of the EFP goes beyond its role as 
an example of such correlators. 

The Emptiness Formation Probability was introduced in the development 
of the determinant representation for correlation functions of the integrable 
models we mentioned before. A good account of this technique can be found 
in [18], where is it shown how the Bethe Ansatz solution can be manipulated 
to produce exact expressions for correlators. A careful analysis of these ex- 
pressions revealed that the simplest correlator one can construct for integrable 
models is precisely the EFP. 

As we mentioned, the practical calculation of correlators for integrable sys- 
tems is still an open challenge. Therefore, it is conceivable that the study of 
the simplest correlator will bring insights helpful to carry on the investigation 
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of other correlators. Probably also for this reason, in recent years a consider- 
able effort has been devoted to the study of the EFP in different systems [H] 
- [30], but a general method for its calculation is still missing. 

We propose the hydrodynamic approach as such a general method for the 
calculation of the EFP, at least to leading order. 

1.4 EFP in different systems 

In this section we introduce the EFP for several simple integrable models. 

Emptiness Formation Probability: Let us first consider a one- dimensional 
quantum system of N particles at zero temperature. The wavefunction of the 
ground state of the system \I/ G (xi,x 2 , . . . , xjv) gives the probability distribu- 
tion \^g\ 2 of having all N particles at given positions Xj, where j = 1, . . . , N. 

The Emptiness Formation Probability P(R) is defined as the probability 
of having no particles with coordinates —R<Xj<R for every j. Formally 
we write this as 




(1.1) 



or following [18] 



P(R)= lim (tf G |e 



f-R p( x ) d: 



X 



(1.2) 
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where p(x) is the particle density operator 

N 

P(x) = ^25(x-Xj). (1.3) 

3=1 

Spin chains and lattice fermions: The EFP can also be defined for spin 
chains. In these systems we are interested in what is known as the "Probability 
of Formation of Ferromagnetic Strings" (PFFS), where we are looking for 
strings of length n in the ground state of the spin chain: 

n 1 z 

^H = <oin^H >' (1.4) 

1=1 

where of is the ^-component of the Pauli matrices on the z-th lattice site. 

The Jordan- Wigner transformation (12.2112. 41) maps a spin- 1/2 chain to a 
one-dimensional lattice gas of spinless fermions. Under this mapping the fer- 
romagnetic string corresponds to a string of empty lattice sites and one can 
write the EFP 

p(R) = ^ n My > ( L5 ) 

where are annihilation and creation operators of spinless fermions on 

the lattice site j. 

Therefore, the PFFS in a spin chain corresponds to the EFP of the cor- 
responding Jordan- Wigner fermion theory. In general, we are going to use 
the language of particles in this work, so we will generically refer to the EFP 
even for spin systems, since all EFP results are valid for the corresponding 
one-dimensional spin systems as well. 
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Random Matrices: As a remark, we point out that the EFP (11.1 h intro- 
duced for a general one- dimensional quantum system is a well-known quantity 
in the context of the spectra of random matrices [36]. Essentially, it is the 
probability of having no eigenvalues in some range of the spectrum. 

Consider, e.g., the joint eigenvalue distribution for the Circular Unitary 
Ensemble (CUE). The CUE is defined as an ensemble of N x N unitary matri- 
ces with the ensemble measure given by the de Haar measure. Diagonalizing 
the matrices and integrating out the unitary rotations, one obtains [61] 

r N 

DU - / J]^ J] \j 6i -j 6k f, (1-6) 

J'=l ^<3<k<N 

where f3 = 2 for the CUE and e l9j , with j = 1, . . . , N, are the eigenvalues of a 
unitary matrix. 

One can then read the joint eigenvalue distribution as 

P N (9 1 ,...,9 N ) = const. je^-e^. (1.7) 

l<j<k<N 

We can now introduce the probability of having no eigenvalues on the arc 
—a < 9 < a as 

. N 



-a,a] 



In the theory of Random Matrices, this quantity is known as £73 (0, a), using the 
notations of |61j . and is clearly equivalent to the EFP once one identifies energy 
eigenvalues and particles. For orthogonal, unitary, and symplectic circular 
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ensembles the joint eigenvalue distribution is given by (11.71) with j3 = 1,2,4 
respectively, where (3 = 2 corresponds to free fermions, while (3 = 1,4 are 
particular cases of the EFP in the Calogero-Sutherland model [61] . 

1.5 Some implementations of One-Dimensional 

systems 

While one-dimensional physics has been an active sector of research in theo- 
retical studies, only recently laboratories around the world have been able to 
build and study devices where the motion of electrons is effectively limited to 
just one dimension and where sufficient control is available to explore different 
regimes and phases. This allows to finally test in experiments the theoretical 
prediction on 1-D models. We present an overview of some systems and some 
experiments where Luttinger liquid behavior has been confirmed and discuss 
whether curvature effects, like the correction to exact spin-charge separation, 
can be observed. 

GaAs/AlGaAs quantum wires: One dimensional quantum wires are re- 
alized by cleaved edge overgrowth on a GaAs/AlGaAs heterostructures. These 
wires are in the ballistic regimes and Luttinger Liquid behavior has been con- 
firmed by several experiments, for instance by measuring deviation from exact 
quantization of the conductance against the theoretical predictions [2]. Re- 
cently, efforts have been devoted in observing signatures of spin-charge sepa- 
ration in these quantum wires [3], but the results are so far non-conclusive. 
The experiment was conducted on two parallel ballistic wires, by measuring 
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the tunneling current: the measurements were in qualitative agreement with 
the theory, but one of the two predicted charge modes was not observed. It 
is conceivable that in the future a similar setup could be implemented to test 
Coulomb drag effects [I], that would carry a clear signs of spectrum curva- 
ture, but attempts in this direction require a level of sophistication beyond 
the technology available today. 

Carbon nanotubes: Single-wall and multi-wall carbon nanotubes are ex- 
cellent examples of effectively one-dimensional systems [51 [6]. Different re- 
alizations of the tubes generated metallic, insulating, semi-metallic and semi- 
conducting wires. Luttinger liquid behavior of metallic carbon nanotubes have 
been demonstrated by measuring tunneling amplitudes into the wires through 
STM experiments j5]. The conductance and differential conductance followed 
power law behaviors as functions of temperatures and bias voltages, and the 
exponents were found to be in good agreement with the theory. For semi- 
metallic wires the description of electrons in terms of Majorana Fermions is 
almost exact in that the spectrum is linear near the conical point with very 
good approximation. In the metallic regimes, the Fermi points move away 
from this point, but we do not expect it to be possible to observe the effect of 
curvature in these systems with the present technology. 

Organic conducting molecules: Some organic molecules, like polyacety- 
lene, can be used as one-dimensional wires. The main limitation in experiments 
with organic molecules is the lack of control in the preparation of the system, 
due to the fact that the molecules have a predefined structure. Therefore, one 
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has to look for the right molecule that could fit a model and there is no a 
generic set-up that can be used to test different regimes. Polyacetylene, for 
instance, is a one-dimensional gapped insulator and therefore it cannot be used 
to test Luttinger liquid behaviors [JJ . Conductance studies on Bechgaard salts 
fibers show the power law behaviors one expects from a Luttinger liquid, but 
the data are hard to quantitatively match with the theory, due to the com- 
plexity of the systems [8]. Evidence for spin-charge separation effects in these 
salts have been gathered by comparing charge and thermal conductivity, but 
for a clear interpretation of the results we would need to know these systems 
better than we currently do. 

One-dimensional metallic chains: For a long time now, one dimensional 
metallic chains like Au atoms on S'z(lll), SrzCuOz and SrCu02S have been 
under investigation, trying to measure evidences of spin-charge separation, 
but results were not conclusive [9]. Recently, a clear signature of spin-charge 
separation has been reported using Angle-Resolved PhotoEmission (ARPES) 
data [10]. They reconstructed the spectral function, which showed two dis- 
tinct peaks corresponding to the spin and charge contributions in quantitative 
agreement with the theoretical prediction. There was, however, a significant 
portion of the spectral function unaccounted by theoretical models, resulting 
in an unexpected broadening of the peaks. It is possible that this effect is due 
to the coupling between spin and charge degrees of freedom arising from the 
spectrum curvature. 
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Cold atomic gases: In recent years, impressing technological developments 
have allowed to cool atomic systems confined in optical traps to unprecedented 
low temperatures. Since the first observation of Bose-Einstein condensation 
[TT] , temperatures have kept dropping, allowing the observation of new physi- 
cal phenomena. An active field of research focuses on atomic gasses effectively 
trapped in one-dimensional geometries. It is believed that many solid state 
physics can be mimicked by these systems and that it would be possible to 
study them at much lower temperatures than their condensed matter counter- 
parts [12] . Most successful have been experiments with bosons, with the recent 
realization of a Tonks- Girardeau gas, i.e. a gas of bosons with such strong 
repulsion that at low enough temperatures and densities exhibit fermionic be- 
havior [13]. It is also expected that a one- dimensional Tonks gas would exhibit 
Luttinger Liquid behavior and some evidence is already available supporting 
this prediction [13]. Direct experiments with fermionic gasses need to decrease 
the temperatures by about two orders of magnitudes compared to the current 
status in order to realize a Luttinger liquid system, but the improvements 
in experimental techniques indicate that this goal is probably not far in the 
future. Indeed, one- dimensional cold atomic gasses are really promising for 
the investigation of curvature effects of the spectrum in the near future, be- 
cause these systems allow for a more direct manipulation of the excitations, 
compared to equivalent solid state systems. 
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1.6 Outlook of the thesis work 



In Chapter [2] we introduce the spin-1/2 anisotropic XY model, an important 
integrable spin chain model in one dimension. In this chapter we analyze the 
model, study its rich phase diagram and derive the fundamental correlators of 
the model. The XY model is arguably the simplest integrable model with a 
non-trivial phase diagram and a perfect candidate for the study of the EFP. 

In Chapter [3] we undertake this study and show that for this system the 
EFP P(n) can be expressed exactly as the determinant of a n x n matrix. This 
matrix belong to a class of very special matrices known as Toeplitz Matrices 
and the asymptotic behavior of their determinant has been studied extensively 
by mathematicians. We use the results of the theory of Toeplitz determinants 
to calculate the asymptotic behavior of the EFP in the different regions of the 
phase-diagram of the model. We find that the EFP decays exponentially in 
most of the phase diagram and only for the isotropic case studied in [23] is 
the decay Gaussian. On the critical lines we observe an additional power-law 
correction to the exponential or Gaussian decay. We employ a bosonization 
approach to interpret the crossover from the Gaussian to the exponential decay 
with universal exponents. These results are original and first appeared in |30j . 

In Chapter H] we introduce the hydrodynamic approach. We first analyze 
the hydrodynamic description of free fermions and then proceed to show the 
technique in its generality and how to obtain the internal energy of integrable 
systems from the Bethe Ansatz. We also show that the linearization of the 
hydrodynamic theory produces the traditional bosonization. 

In Chapter [5] we show how to derive the leading asymptotic behavior of the 
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EFP for some Galilean invariant systems. This calculation is performed at a 
semi-classical level, where the EFP is calculated as the probability of a rare 
fluctuation ( "instanton") . We demonstrate that bosonization is insufficient to 
quantify the EFP. Using the hydrodynamic approach we calculate the EFP to 
leading order for Free Fermions and for Calogero-Sutherland particles. 

In Chapter [6] we analyze the prediction of spin-charge separation in one 
dimension. First, we show the limitations of the bosonization approach in 
dealing with the correction to exact spin-charge separation. Then we develop 
a hydrodynamic description for fermions with contact repulsion. While our 
ultimate goal is to address the problem of evaluating the corrections to exact 
spin-charge separation, this is still a work in progress and we can only show 
how to implicitly construct the hydrodynamic description of this model from 
its Bethe Ansatz solution. 

In Chapter [7] we address a problem very different from the rest of the thesis. 
It is a single-particle problem in two dimensions. We study the behavior of 
a scalar particle in a medium filled with point-like magnetic fluxes (vortices) 
pinned on the sites of a square lattice. By assuming the strength of the fluxes 
to be all equal to half of the quantum flux unit, we are able to construct a 
wavefunction of the system, showing that the spectrum of such a system is 
discrete. Moreover, we are able to show that a zero-energy particle entering 
such a medium would decay exponentially. This means that an array of vortices 
could be used to trap a particle and this bound state would have a topological 
nature. As far as we know, this is the first time that the existence of such a 
bound state is suggested. 

Finally, in Chapter M we will discuss the results of this thesis, the problems 
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that remain open and the directions for further research. 

We include appendices to supplement the analysis of the main text. 

Appendix |X] is a recapitulation of known results on the EFP calculated for 
various integrable models. 

Appendix [B] contains a brief review of the results of the theory of Toeplitz 
determinants, which is extensively used in the evaluation of the asymptotic 
behavior of the correlation functions of the XY model, including the EFP we 
study in Chapter [3j 

Appendix O is a very brief introduction to the Bethe Ansatz technique 
where we gather some results we need as an input for the hydrodynamic de- 
scription of integrable models. 

In Appendix [D] we construct the conserved quantities of gradient-less hy- 
drodynamic theories. The integrability of Galilean invariant hydrodynamic 
models without gradient correction is a little-known fact discovered in the 
1980's. Here we explicitly construct the integrals of motion for the first time 
and we attempt to clarify the origin of this integrability. 
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Chapter 2 

The Spin-1/2 Anisotropic XY 
Model 

The One-Dimensional Spin-1/2 Anisotropic XY spin chain is arguably one 
of the simplest non-trivial quantum integrable models. The reason for this 
simplicity lies in the fact that it can be reduced to a free fermions model. 

The XY model describes a one dimensional lattice system, where each lat- 
tice site is occupied by a spin-1/2 quantum degree of freedom interacting with 
its neighbors. The allowed interaction involves only the X and Y components 
of the spins. In addition, we consider the presence of an external transverse 
magnetic field interacting with the Z component of the spins. 

Using what is known as a u Jordan- Wigner Transformation" the spin de- 
grees of freedom can be mapped into spin-less fermions, so that the model 
describes lattice fermions with a quadratic Hamiltonian. The Hamiltonian be- 
ing quadratic, a "Bogoliubov Rotation" defines "Bogoliubov Quasi- Particles" 
in terms of which the model reduces to lattice free fermions. 
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The apparent simplicity of the model does not mean that it is trivial, as 
these quasi-particles are non-local in terms of the original degrees of freedom. 
Therefore, every quantity that one wants to calculate for the original model 
has a non-trivial expression in terms of the free fermions. 

A great simplification comes from the fact that for the XY model these 
non-trivial expressions can often be expressed as determinants of matrices with 
remarkable symmetries. These matrices are known at "Toeplitz Matrices" and 
a rich mathematical literature has been devoted to the study of the asymptotic 
behaviors of their determinant. We will make extensive use of these results on 
Toeplitz determinants. 

The anisotropic XY spin- 1/2 chain in a transverse magnetic field is defined 
by the Hamiltonian 

-hJ2«i, (2-1) 

i=l 

where erf, with a = x,y, z, are the Pauli matrices which describe spin operators 
on the z-th lattice site of the spin chain and, for definiteness, we require periodic 
boundary conditions: erf = crf +N (N » 1). 

We are interested in this model because of its rich phase-diagram. The 
model depends on two parameters, the external magnetic field h and the 
anisotropy parameter 7 controlling the relative strength of the interaction be- 
tween the X and Y components of the spins. By varying these parameters the 
system crosses several quantum phase transitions (QPT). The rich structure 
of its phase diagram is in contrast with the fair simplicity of the model and 
makes it an excellent candidate for examination in connection with QPT. 



A — 1 



1+7 



<J : a. 



i+1 



+ 



7 



v v 
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This is one of the reasons for which the XY model is one of the most studied 
in the growing field of quantum computing and quantum information. 

In section 12.11 we introduce the model, and show its different formulations 
in terms of spins and spinless fermions, and we analyze its phase-diagram. 
In section 12.21 we derive the fundamental correlators for the model. In this 
section we will also give an overview of the results of McCoy and co-authors 
[3T] to better understand the structure of this model in preparation for the 
next chapter where we will calculate the EFP for the XY model [2U [30] . 

The XY model was first solved by Lieb, Schultz and Mattis in [32] without 
the external magnetic field. McCoy and co-authors [31] were the first to de- 
velop the theory of Toeplitz determinants in connection with the study of the 
XY spin chain in the external field and to exploit these structures to derive 
the fundamental correlators of the model. 



2.1 The model and its phases 

The XY model described in (12. ip has been solved in [32] in the case of zero 
magnetic field and in [31] in the presence of a magnetic field. We follow the 
standard prescription [32] and reformulate the Hamiltonian (12. ip in terms of 
spinless fermions ipi by means of a Jordan- Wigner transformation: 



j*Z k< A^ = $ H (2^ k - l) , (2.2) 

k<j 

e -i*E*<^ = ft {2^1 - l) Vi, (2-3) 



k<j 



a] = 2^-1, (2-4) 
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where, as usual, 



a ± = . (2.5) 



The Hamiltonian in terms of these spinless fermions becomes: 

N 

H = Y1 ty^+i + ^+1^ + 1 + i - 2h ^M) ( 2 - 6 ) 

1=1 

and in Fourier space it reads (ipj = Ylq^q^)'- 

H = \2 (cosg - ft) -0g^ 9 + 17 sin ? V^-g - *7 sin <? ^-<?^<z • ( 2 -7) 
q 

We can now diagonalize this Hamiltonian by means of a Bogoliubov trans- 
formation 



X q = cos^ ^ + i sin-^ i/>t q , (2.8) 



which mixes the Fourier components with "rotation angle" d q defined by 



£ 9 



■(cos g — h + i7 sin q). (2.9) 



In terms of these new Bogoliubov quasi-particles Xq the Hamiltonian (12.71) 
has diagonal form 

# = X>xJx* ( 2 - 10 ) 



q 

with the quasiparticle spectrum 



e q = y (cos q — ft) 2 + 7 2 sin 2 g. (2-11) 
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Figure 2.1: Phase diagram of the XY Model (only the part 7 > is shown). 
The theory is critical for h — ±1 (fi±) and for 7 = and \h\ < 1 (fio)- The line 
Ti represents the Ising Model in transverse field. On the line Te the ground 
state of the theory is a product of single spin states. 
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A system is said to be "critical" when its spectrum is gapless, i.e. when 
one can excite particles from the Fermi Sea without spending energy. When a 
system becomes critical, it undergoes a "Quantum Phase Transition (QPT)". 
QPTs are zero-temperature analogs of traditional phase transitions. QPTs 
are characterized by singularities in thermodynamic quantities and by corre- 
lators having a characteristic algebraic behavior. The effective theory is scale 
invariant and in one dimension can be described through Bosonization. QPTs 
are a very active field of research, especially since experiments can reach low 
enough temperatures where their signatures are observable. 

We recognize from (12.111) that the theory is critical, i.e. gapless, for h = ±1 
or for 7 = and \h\ < 1. 

In Fig. I2.1l we show the phase diagram of the XY model, which has obvious 
symmetries 7 — > —7 and h — > — /o The phase diagram has both critical and 
non-critical regimes. Three critical lines f2 (Isotropic XY model: 7 = 0, 
\h\ < 1) and Q± (critical magnetic field: h = ±1) divide the phase diagram 
into three non-critical domains, E_, S , and S + (h < —1, —l<h<l, and 
h > 1 respectively). Fig. 12.11 also shows the line 7 = 1 (Tj) corresponding 
to the Ising model in transverse magnetic field and the line 7 2 + h 2 = 1 (r^) 
on which the wavefunction of the ground state is factorized into a product of 
single spin states 



1 In the next chapter we will analyze the EFP for the XY and we will see that this 
correlator breaks the latter symmetry. 
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2.2 The correlators of the model 



In this section we review the derivation of the fundamental correlators for the 
XY model at zero temperature following McCoy and co-authors [3T] . 

The ground state |0) of the model in terms of the Bogoliubov quasi-particles 
in f l2.10p is defined as 

X g |0> = Vg (2.12) 

i.e. it is the conventional ground state of free fermions. The correlators for 
this theory are easily found to be 

(0\x\Xk\0) = 5 k>q , (2.13) 

(0|X 3 XI|0) = 0, (2.14) 

(0\x q Xk\0) = 0, (2.15) 

<0|4x!l°> = 0. (2.16) 

This vacuum is the ground state for the XY model, but it is not so simple 
when expressed in terms of physical particles. The Hamiltonian (I2.7P contains 
superconducting-like terms, so its ground state is non-trivial. One can invert 
the Bogoliubov transformation ( 12.81) 



^ = cos^Xg-isin^X-g ( 2 - 17 ) 



to calculate the fundamental correlators in terms of physical fermions: 



(0|^ fc |0) = sin 2 ^5 M = 1 C 2 ° S *H , g , (2.1* 
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<0|^!|0> = cos 2 !^ = 1 + C 2 ° S ^H , g , (2.19) 

$ $ sin 

(0|^fc|0) = -i cos -j- sin -j- 5 kjq = -i S k>g , (2.20) 

(o|^Il°) = icos y^y^ = i^^- (2-21) 



Now, the two-point fermionic correlators are easy to obtain by Fourier 
transform. In the thermodynamic limit they read [32j [31] 



F jk ee mMk\0) = -m$4\V) = J*^^^-*\ (2-22) 
G jk - (0|^|0) = ^|i±f^e^-). (2.23) 

These correlators will be fundamental in our calculation of the EFP in the 
next chapter. 

To calculate the correlation functions for the original spin chain model 

PL = (0K<J0> u = x,y,z, (2.24) 

we need more work. We follow [32] and express these correlators in terms of 
spin lowering and raising operators (12.51) : 

Pf m = (0|(tf + fff) K + a")|0), (2.25) 
PL = (0|K"^r) K-(T-)\0), (2.26) 
p\ m = (0 | (2a+ ff f - 1) (2<x> m - 1) | 0). (2.27) 

and use (I2.2H2. 41) to write them using spinless fermions operators. 
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For instance, let us consider pf m : 

Pirn = (OlK + ^r) (^m + ^m)|0> 

m— 1 

= (o| U\ + tin) [ [ Ml* - i) (Vi + <M 1°) 



i=l 
m—1 



(oi(^?-^) n (^m- 1 ) w+^)io) 
t=i+i 

m— 1 

(oi (W - ^) n + ^) - ^) ml + ^) ( 2 - 2 ^ 



i=Z+l 



where we have used two identities 



a+ = e ^<i ^ fc = e - i7rEfe <^ ^ fc ipj (2.29) 



and 



e i^M = _ ! = + ^ j _ ^ = _ _ ^ + ^ (2 30) 

Now we define the operators 

Ai =4 + ^i (2-31) 
5, =4-^ (2-32) 

which allow us to write the correlators (I2.24p as 

pf m = (U\B l A l+1 B l+1 ...A m _ l B m _ x A m \ti) 

Pin = (-l) m - 1 (0\A l B l+1 A l+1 ...B m _ 1 A m _ 1 B m \0) 
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pf m = (0|A,B,AA|0). (2.33) 

We can use Wick's Theorem to expand these expectation values in terms 
of two point correlation functions. By noticing that 

(0|^An|0) = (0|^£y0> = (2.34) 

we write pf m as 

p z lm = (0|A^|0)(0|A mJ B m |0) -(oi^s^oxoia^^io) 

= H 2 (0)- H(m-l)H(l-m) (2.35) 

where 

H(m -l) = (0|AA m |0) = \ TP- e^ m ~ l \ (2.36) 

2 J lix 

The other two correlators in (I2.33H involve more terms. It can be shown [3T1 
[32] that the Wick's expansion can be expressed as a Pfafhan of a matrix with 
elements given by expectation values of each combination of two operators. 

The Pfaffian of a matrix M is defined as [33J 

Pf(M) = ^(-1) P M P1P2 M P3P4 . . . M P2n _ lP2n , (2.37) 
p 

where P = {pi,P2, ■ ■ ■ ,P2n} is a permutation of {1, 2, . . . , 2n}, the sum is 
performed over all possible permutations, and (— l) p is the parity of the per- 
mutation. 
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By using one of the fundamental properties of the Pfaffian: 



Pf(M) = v /det(M) (2.38) 
we can write the spin correlators (12.331) as m — I x m — I matrix determinants: 



PI 



Pirn 



(2.39) 
(2.40) 



with matrix elements given by (12.361) . 

Matrices like (12. 39112. 40p are very special. Their entries depend only on 
the difference between the row and column index, so that the same elements 
appear on each diagonal. Therefore they look like: 





H(-l) 


H{-2) H{- 


-3) 




H(- 




H(0) 


H(-l) H(- 


-2) 




H(l 


Pfm = P X (N) = 


H(l) 


H(0) H(- 


-1) 




H(2 




H(N -2 


) H(N-3) H(N 


-4) 




H(- 




H(l) 


H(0) H(-l) 




H(2 


-N) 




H{2) 


H(l) H(0) 




H(3 


-N) 


PL = P V (N) = 


2f(3) 


H{2) H(l) 




H(A 


-N) 



(2.4i; 



(2.42) 



H(N) H(N-l) H(N-2) ... H(l) 
where N = m — I. 

Matrices like (12.39112.401) are known as "Toeplitz Matrices" and a vast 
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mathematical literature has been devoted to the study of the asymptotic 
behavior of their determinants ( "Toeplitz Determinants") . McCoy and co- 
authors [31] were among the first to develop the theory of Toeplitz determi- 
nants in connection to physical systems. They showed that the fundamental 
correlators for the XY model can be calculated in terms of Toeplitz Determi- 
nants. In the next chapter we are going to calculate the Emptiness Formation 
Probability in the XY model, a non trivial correlator, and we are going to 
show that it also can be expressed as a Toeplitz determinant. 

We will defer an appropriate discussion on the theory of Toeplitz matrices 
to the next chapter, where we are going to use it to calculate the EFP. Here, 
we are just going to recap the main results of [31] on the asymptotic behavior 
of the fundamental spin correlators (12. 24ft . 

At zero temperature, the asymptotic evaluation of the Toeplitz determinant 
in the different regions of the phase diagram for p x (N) gives [31] 



p^Nf- 00 



(-irT^[7 2 (l-^)] 1/4 |*|<1,7>0 

^l)N^_ e i/i 2 i/i2 A -3^ N yi/i [ft] = 1, 7 > o 

(-l)"/i(*>7)$ W>1,7>0 

10 7 = 



(2.43) 



where A = 1.282. . . is the Glaisher's constant, fi(h,j) is some function (see 
[31]) and 

>,^vgHI. (2.44) 
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For p y (N) it was found 



r -(-i) N h(h n )N-*\-™ 





p y(N) N ~°° { 



\h\<^l^f ll >0 



■l) N h(h, 7 )N-i fe) 



N 



\h\ = v/r^, 7 > o 

Vl+7 y VT^<|/1|<1, T >0 

(_ 1 )AT7(l+7l e l/4 2 l/12 A -3( 7Ar )-9/4 \h\ — 1, J > 



-1) N f i (h,'j)N~ 3/2 \ 



N 



I 



Finally, for p z (N) we have 



\h\ >1, 7 >0 
7 = 0. 



(2.45) 



p (N) ~ 



m; 



/s(/i,7)^ 



-2 f 1^2 
1+7 



2V 



m 2 - ^A~ 2 ^ 2 
m 2 - (nN)- 2 



1 - iiV-A 



-2 AT 



m z 
i 

I 4 



Sm(7VcOS lh) 



\h\< ^1 "7 2 , 7>0 
= Vl-7 2 , 7 > 



1 -7 2 < \h\ < 1, 7 > 
|/i| = l,7>0 
|/ i |>l, 7 >0 
|/i|< 1,7 = 
|/i|> 1,7 = 



(2.46) 



where m 2 is the magnetization: 



h — cos q 



\J {h — cos q) 2 + 7 2 sin 2 g 



dg. 



(2.47) 
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Chapter 3 



The EFP for the XY Model 

In the previous chapter we introduced the One Dimensional Spin-1/2 Anisotropic 
XY spin chain in a transverse magnetic field. We studied its phase diagram 
and we calculated the fundamental correlators of this model. Now we turn our 
attention to a non-trivial correlator know as Emptiness Formation Probability 
(EFP) which we are going to study for the XY model. 

We introduced the EFP in Chapter 1, where we discussed its significance 
and importance in the theory of integrable models and in the general problem 
of calculating correlators in one-dimensional theories. 

The XY model is a very interesting model for the study of the EFP for 
several reasons. As we saw in the previous chapter, the XY spin chain is char- 
acterized by a very interesting phase diagram: as we vary its two parameters 
(the anisotropy and the magnetic field), we move through critical and gapped 
regions. Previous studies of the EFP focused only on critical phases of various 
systems and the XY model offers an opportunity to follow the behavior of the 
EFP across a phase transition. 
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A second, more technical reason, lies in the relative simplicity of the model. 
In the previous chapter we showed that the fundamental correlators of the 
theory can be exactly expressed as determinants of a very special class of 
matrices known as "Toeplitz Matrices". We show that the same property 
holds for the EFP. This is of great advantage, since, despite being considered 
the simplest correlator for integrable models, the EFP in general does not have 
a simple expression. Using the theory of "Toeplitz Determinants", we are able 
to calculate the asymptotic behavior of the EFP in the various regions of the 
phase diagram of the XY model. 

Most of these results appeared first in [30]. The Toeplitz determinant 
approach was also used in Ref. [24J for the EFP in the case of the Isotropic XY 
model (Eq. ( 12. ID with 7 = 0). In this latter work it was shown that the EFP 
decays in a Gaussian way for the critical theory (7 = 0, — 1 < h < 1). This 
case corresponds to one of the two critical lines in the 7-/1 phase diagram of 
the model ( 12. ID discussed in the previous chapter. The other line is the critical 
magnetization line(s) (h = ±1). In the rest of the two-dimensional 7-/1 phase 
diagram, the model is non-critical. 

For the XXZ spin chain in zero magnetic field, the EFP was found to have 
a Gaussian decay P(n) ~ e~ an2 as n — > 00 in the critical regime at zero 
temperature and exponential e _/3n at finite temperature (|28|.[29j). 

A qualitative argument in favor of Gaussian decay was given in Ref. [23] 
within a field theory approach. It was argued there that the asymptotics 
of the EFP are defined by the action of an optimal fluctuation (instanton) 
corresponding to the EFP. In the critical model, this fluctuation will have 
a form of an "n x n" droplet in space-time with the area A ~ n 2 and the 
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corresponding action S ~ an 2 which gives the decay Pin) ~ e~ an2 . Similarly, 
at finite temperature the droplet becomes rectangular (one dimension n is 
replaced by an inverse temperature T _1 ) and the action cost is proportional 
to n, giving P(n) ~ e~@ n . This argument is based on the criticality of the 
theorjo and it is interesting to consider whether it could be extended to a 
non-critical theory A naive extension of the argument would give the optimal 
fluctuation with space-time dimensions n x £ where £ is a typical correlation 
length (in time) of the theory. This would result in P(n) ~ e~ /3n for non- 
critical theories, similarly to the case of finite temperature in critical regime. 
The rate of decay (3 would be proportional to the correlation length of the 
theory. 

In this chapter we examine the relation between the asymptotic behavior 
of the EFP and criticality using the example of the XY model. Using Toeplitz 
determinant techniques, we obtain that the EFP is asymptotically exponential 
in most of the phase diagram according to the naive expectations and that it is 
Gaussian only at 7 = in agreement with previous works on XXZ spin chains 
and with Ref. [21]. However, on the critical lines h = ±1, in addition to the 
exponential decay, a pre-exponential power-law factor arises, with a universal 
exponent. The power-law prefactor is present in the isotropic CclSG clS well, but 
with a different exponent. Using a bosonization approach, we will interpret 
the transition from Gaussian to exponential decay. 

The chapter is organized in the following way: in Section [3] we explain how 
one can express the EFP as the determinant of a Toeplitz matrix and review 



^^Morc precisely, on the assumption that temporal and spatial dimensions of an instanton 
scale similarly. 
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our results so that readers who are not interested in derivations can skip the 
next sections. In Section 13.21 we analyze the exponential decay of the EFP for 
the non-critical and critical phases of the anisotropic XY Model. In Section 
I3.3l we derive in detail the asymptotic behaviors, including the pre-exponential 
factors, of both non-critical and critical parts of the phase diagram. In Section 
13.41 we study a special line of the phase diagram on which the ground state 
is known exactly and compare the explicit results one can obtain using the 
exact ground state with the asymptotes of the EFP we derived in the previous 
sections. In Section 13.51 we report on the already known results for the EFP 
of the isotropic XY model [21] . In Section 13.61 we make contact with Ref . 
[23] using a bosonization approach to discuss the crossover as a function of n 
from the Gaussian to the exponential behavior of EFP for the case of small 
anisotropy 7. The following section gives some mathematical details on the 
calculation of the stationary action in the bosonization approach of Section 
13.61 and can be skipped by the reader non interested in the mathematical 
technique. Section [3781 presents the analysis of the finite temperature behavior 
of the EFP, which gives an expected exponential decay. Finally, Section I3~9l will 
summarize our results. For the reader's convenience we have collected some 
results on asymptotic behavior of Toeplitz determinants which are extensively 
used in the rest of the paper in the Appendix OH 

3.1 EFP as a determinant of a Toeplitz Matrix 

We introduced the XY model in the previous chapter. In this section we con- 
sider the correlator measuring the "Probability of Formation of Ferromagnetic 
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Strings" (PFFS) 

n 

at zero temperature (the non-zero temperature case is deferred to section 



i=l 



13. 8p . This correlator measures the probability that n consecutive spins will 
be aligned downwards in the ground state of the system. 

In terms of the spinless fermions defined in the previous chapter (12.71) . by 
direct substitution one can express the PFFS (13.11) as the expectation value 
over the spinless fermions ground state [24] 



PH = (0|n^J|0). (3.2) 



This expression projects the ground state on a configuration without particles 
on a string of length n and hence gives the meaning to the name "Emptiness 
Formation Probability". 

Let us now introduce the 2n x 2n skew-symmetric matrix M of correlation 
functions 

f-iF G\ 

M = , (3.3) 

V-G iFJ 

where F and G are n x n matrices with matrix elements given by Fjj. and 
Gjk from (12.22112.231) respectively. Then, using Wick's theorem on the r.h.s of 
(13. 2p . we obtain the EFP as the Pfaffian of the matrix M 



P(n) = Pf(M), (3.4) 
where the Pfaffian was introduced in (12.371) . Using one of the properties of the 
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Pfaffian we have 



P(n) = Pf (M) = v / det(M). (3.5) 
We can perform a unitary transformation 

/ S n \ ! (I -I\ 

M' = UMU f = , U = -= , (3.6) 

V-s n t Oj V2\l I J 

where I is a unit n x n matrix and S n = G + iF and S n ^ = G — iF. This 
allows us to calculate the determinant of M as 

det(M) = det(M') = det(S„) • det(S n f ) = |det(S n )| 2 . (3.7) 

The matrix S n is a n x n Toeplitz matrix (i.e. its matrix elements depend only 
on the difference of row and column indices [37] like the matrices defined in 
the previous chapter f l2.39U2.40f) ) . The generating function a(q) of a Toeplitz 
matrix is defined by 



>n)jk 



2 " ^ <T(q)<j*V-Q (3.8) 



and in our case can be found from (I2.22ll2.23j) as 

a(q) = 1(1 + e^) = \ + ^j^A+j^ . (3 . 9 ) 
2 2 2 a/ (cos q — h) 2 + r y 2 sin q 

Thus, the problem of the calculation of the EFP 

P(n) = |det(S n )|, (3.10) 
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is reduced (exactly) to the calculation of the determinant of the nxn Toeplitz 
matrix S n defined by the generating function (13. 8113. 9p . The representation 
(I3.10P is exact and valid for any n. In our study we are interested in finding 
the asymptotic behavior of (I3.10p at large n — > oo. |^| 

Most of these results are derived using known theorems on the asymp- 
totic behavior of Toeplitz determinants. We have collected these theorems in 
Appendix [B] In the following sections we apply them to extract the corre- 
sponding asymptotes of P(n) at n — > oo in the different regions of the phase 
diagram. Two major distinctions have to be made in this process. For the 
critical isotropic (7 = 0) XY model, one applies what is known as Widom's 
Theorem and finds a Gaussian behavior with a power law prefactor [23] . In the 
rest of the phase diagram, we apply different formulations of what is known in 
general as the Fisher-Hartwig conjecture, which always leads to an exponential 
asymptotic behavior. As expected, we find a pure exponential decay for the 
EFP in the non-critical regions. 

For h > 1, the exponential decay is modulated by an additional oscillatory 
behavior. 

At the critical magnetizations h = ±1, we discover an exponential decay 
with a power law pre-factor. Moreover, by extending the existing theorems 
on Toeplitz determinants beyond their range of applicability, for h — ±1 we 
obtain the first order corrections to the asymptotics as a faster decaying power 

law with the same exponential factor. For h — 1, the first order correction is 

2 The reader might notice that our generating function (|3.9p is almost the same as the one 
analyzed by Barouch et al. in [31] (<J\ 13] (q) = cos <i- h +h si "<? ). The only difference is 

L J w (cos q— /i) 2 +7 2 sin 2 q 

the shift by the unity in our expression. This difference changes dramatically the analytical 
structure of the generating function, in particular, its winding number around the origin, 
and requires a new analysis of the generated Toeplitz determinants. 



36 



h 



o 







[ 1 + A cos(Tin) ] e 


o 

1 Q 

4- 




Bn 

, n V16 [ 1 + A n 1/2 cos(7in) ] e 


I 

2 

ri /4 e 


— \— 

Crossover a 


'V 
I 


* ► e 






s — 


-pn 

e 



Figure 3.1: Asymptotic behavior of the EFP in the different regions of the 
phase diagram of the XY Model (only the part 7 > is shown). The theory is 
critical for h — ±1 (fl±) and for 7 = and \h\ < 1 (fin)- The line Tj represents 
the Ising Model in transverse field. On the line Te the ground state of the 
theory is a product of single spin states. 

also oscillating, which indicates that the EFP has an oscillatory behavior for 
h>l. 

The reader who is not interested in the mathematical details of our deriva- 
tions can find the results in Figure [37T1 and in Table [3781 and skip the following 
sections to go directly to Sec. 13. 6} where we analyze the crossover between the 
Gaussian behavior at 7 = and the asymptotic exponential decay at finite 7 
using a bosonization approach. 
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Figure 3.2: Plot of the 
_ _- — — absolute value and argu- 
ment of the generating 
function (13. 9p for 7 = 1.5 
at different values of h. 
From (a) to (e) h = — 1.1, 
—1, 0.5, 1, 1.1, respec- 
tively. 



3.2 Singularities of a(q) and exponential be- 



havior of the EFP 



To derive the asymptotic behavior of the EFP we rely on the theorems on 
determinants of Toeplitz matrices. These theorems depend greatly on the 
analytical structure of the generating function (13.91) . especially on its zeros 
and singularities. 

Setting 7 = in (13.91) . we see that for the Isotropic XY model the gener- 
ating function has only a limited support within its period [0, 2tt]. This case 
is covered by what is known as Widom's Theorem and will be considered in 
Section 13.51 

In the remaining parts of the phase-diagram the generating function has 
only pointwise singularities (zeros) as it is shown in Fig. 13.21 These cases 
are treated under a general (not yet completely proven) conjecture known as 



38 



the Fisher-Hartwig conjecture (FH), which prescribes the leading asymptotic 
behavior of the Toeplitz determinant to be exponential in n: 




(3.11) 



While the pre-exponential factors depend upon the particulars of the sin- 
gularities of the generating function, the exponential decay rate is given in the 
whole phase diagram (7 ^ 0) according to FH as 



The integral in (13.1 2p is convergent for all h and all 7 7^ and /3(h,j) is a 
continuous function of its parameters. 

In Fig. I3.3[ (3(h, 7) is plotted as a function of h at several values of 7. One 
can see that fl(h, 7) is continuous but has weak (logarithmic) singularities 
at h = ±1. This is one of the effects of the criticality of the model on the 
asymptotic behavior of EFP. 

These weak singularities are also a manifestation of the rich analytical 
structure underlying /3(h, 7) and the generating function (13.91) . To understand 
these structures, we interpret the periodic generating function (13. 9p as the 
restriction to the unit circle (z = e l9 ) of the complex function 




(cos q — h) 2 + 7 2 sin 2 q 



cos q — h 



.(3.12) 




(3.13) 
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Figure 3.3: Plot of the decay rate (3 as a function of the parameters 7 and h. 
The function diverges for 7 = and is continuous for h — ±1 (although it has 
weak singularities at h = ±1). 



where 



P2O) : 



1+7 
" 2z 
1+7 



(z - zx)(z - z 2 ), 



2z 



{z x z - l)(z 2 z - 1) 



(3.14) 
(3.15) 



with 



~2 



h-y/h?+T 2 

1+7 



h + y//i 2 + 7 2 - 1 
1 + 7 



(3.16) 
(3.17) 



The integral in (13.121) can be regarded as a contour integral over the unit 



10 




Figure 3.4: The integral in (13.121) is performed over the unit circle C . The 
analytical structure of the integrand allows for a deformation of the contour 
of integration into T, which encloses a logarithmic branching line, different in 
the various regions of the phase-diagram (in (d), T encloses also a simple pole 
at the origin). The roots z\ and Z2 were defined in ( 13 . 16f) and (13. 17ft . 



circle of the function (13.131) . We can deform the contour of integration taking 
into account the complex structure of the integrand in the various regions (see 
Fig. I3.4p and express (13.121) as a simpler integral on the real axis (after partial 
integration and some algebra). 
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3.2.1 The non-critical regions (£± and So) 
3.2.1.1 £_ (h < -1) 



For h < — 1, the analytical structure of the integrand of (13.121) is shown in 
Fig. 13.4b . We re-write the decay rate A3. 12|) in this region as 



In 



^h 2 + 7 2 - 1 - ft 
7 + 1 



A(/i, 7 )-A(/i, 7 ), 



(3.18) 



where 



A(fc, 7 ) = In 
A(/i )7 ) = 



1 



1 



h 1 -7 



2 V" WVl + 7 
1 dx 1 

|AT| 2^v^T^ 2 )(x 2 -K 2 ) 







x — a 












x + a 



(3.19) 
(3.20) 



with 



K = 



r 



7 



V^ 2 + 7 2 



1+7 
T-7 



h-1 



(3.21) 
(3.22) 



This decomposition of /3(h, / -f) is especially useful in analyzing the transitions 
between non-critical and critical regimes. In fact, we will see that the func- 
tions A(h, 7) and A(h, 7) defined above are universal across the phase diagram 
(hence the need for the seemingly redundant absolute values in our defini- 
tions). The absolute value in the logarithm of the integrand is relevant for 
7 > 1, since its argument changes sign and vanishes within the interval of 
integration (a > K for 7 > 1). 
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3.2.1.2 S (\h\ < 1) 

As before, we can express the contour integral defining fl(h,j) as a standard 
integral on the real axis. For \h\ < 1 and h 2 + 7 2 > 1, the structure of the 
integrand is depicted in Fig. 13.4b and the decay rate is simply 



/3(/i, T ) = -A(/i, T )-A(/i, T ), (3.23) 

where A(h,/y) and A(/z,7) have already been defined in (13.19113.201) . 

For h 2 + 7 2 < 1, the structure is quite different (see Fig. I3.4d ). In this 
region the expression for /3(h, 7) in terms of a real axis integral is complicated. 
We can write the result as: 



/?(/*, 7) = -A(/i,7)-v / 2(l + 7)(l-^ 2 -7 2 )x 



^ dx 
— ( 







x 
lh 



V^-h 2 - 



A{x) + B(x) + C{x) 



-\h\ X 



(1+7) 2 



VI - x 2 sj y/q(x) +p{x) i 

A(x) - B(x) - C(x) 
yi - x 2 \J \/q(x) -p(x) i 



(3. 



where A(/i, 7) was defined in (13.191) and 



;) = (I+7) 



r(x) 



t(x) 



S (X) 



A(x 

B(x) = 2(1 + 7) 
C (x) = ( 7 2 -i + 2 7 /i 2 ) 



t(x)v/g(x) 



(3.25) 
(3.26) 
(3.27) 
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with 



p(x) 
q(x) 



r(x) 
s(x) 

t(x) 



(7 + l) 3 - (3 7 + l)h 2 + (7 - 1)(1 - h 2 - 7 V, (3-28) 
(7 + l) 6 - 2(37 4 + IO7 3 + 12 7 2 + 67 + l)h 2 + (97 2 + 67 + l)h 4 
+2[ 7 4 + 2 7 3 - 2 7 - 1 + (57 2 + 2 7 + l)h 2 }(l - h 2 - -? 2 )x 2 
+ ( 7 -l) 2 (l-/i 2 -7 2 )V, (3.29) 
(7 + l) 2 - (2 7 + l)h 2 + (l-h 2 - 7 V, (3.30) 
(7 + l) 4 - (3 7 4 + 8 7 3 + 9 7 2 + 67 + 2)h 2 + (5 7 2 + 2 7 + l)h 4 
+ [(7 + l) 2 - (3 7 2 + 67 - l)h 2 ](l -h 2 - 7 2 )x 2 , (3.31) 
[(7 + 1 + hf + (l-h 2 - 7 2 )x 2 ] [(7 + 1 - h) 2 + (l-h 2 - 7 V]. 

(3.32) 



3.2.1.3 £+ (h > 1) 



A calculation similar to the previous ones (see Fig. I3.4F ) gives the expression 
for the decay factor for h > 1: 



= -In 



v//i 2 + 7 2 - 1 + h 
7 + 1 



-A(h n )-A(h n ), 



(3.33) 



where A(/i, 7) and A(/i, 7) were introduced in (13 . 1 9[) and (13.201) . 

One important difference exists in this region: as will be discussed in length 
later in Section 13.3.1.31 in S + there are two equivalent representations of the 
generating function. This ambiguity reflects on the value of f3, in that the 
choice of the representation for the generating function determines the branch 
cuts in Fig. 13.41 We will see that we have to use both values of f3, which differ 
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only by an imaginary constant: 



(3' = (3 + in 



(3.34) 



and this will add an oscillatory behavior to the EFP. 
3.2.2 The critical lines (Q±) 

We can calculate the decay factor (3 at h — 1 (f2+) from a limiting procedure 
on (I3.23P or (I3.33p . At h = 1, only A(h,j) is nonvanishing, thus guaranteeing 
the continuity of (3 across the critical line. From an appropriate limit of (I3.20p . 
we calculate the decay rate for h = 1 as 



0(1,7) 



1 dx 



In 



1 — 7X 



1 + 7X 



In 



3' 



1-7 
1 + 7 



(3.35) 



For 7 < 1, we can expand the logarithm in series and perform the integral: 



/3(1, 7 <1) = -p 



j oo 



1 



2n+l 



^^r(n + l/2) (2n + l) 2 



In 



2V Vl+7 



. (3.36) 



As discussed before in reference to (13.331) . the definition of (3(1, 7) is not unique 
and, as in (13.341) . will generate again an oscillatory behavior for the EFP (see 
later in Sec. 13.3.2.11) . 

The value of (3 at h — — 1 can also be obtained from a limiting procedure 
on (gj20j) 



0(-l,7) 



1 dx 1 

2^"7r^2 



In 



1 — 72 
1 + 7X 



-In 



1-7 
1+7 



(3.37) 
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Again, for 7 < 1, we can expand the logarithm in series to calculate the 
integral: 



^ ^ < ^ = ~ ^ £ r(n + l/2)(2n + l )a " ^ [2 V + V ~ J 



(3.38) 



As can be seen in Fig. 13.31 the decay factor (3 is continuous across the critical 
lines, but has a discontinuity of its derivative. As (3 approaches the critical 
lines, it actually shows a non-analytical behavior leading to a logarithmic 
singularity: 

= ±1 + e, 7) =/3(±l, 7 ) + -eln|e|. (3.39) 

7T 

The derivative dj3/dh diverges logarithmically as h — > ±1. 

Moreover, one can easily notice from the difference between expression 
(I3.23P and f!3.33j) that even the finite part of the derivative of fl(h, 7) by h is 
different if one approaches the critical line h = 1 from above or below, due 
to the appearance of the additional term in (13.331) . The same holds across 
the critical line h = — 1, due to the presence of the first term in (13.181) . which 
doesn't appear in 1 13.231) . 



3.3 The pre-exponential factors 

For 7 7^ 0, the leading behavior of the EFP is always exponential. However, 
the singularities of the generating function are different in different regions 
of the phase diagram and we must therefore use different forms of the Fisher- 
Hartwig conjecture in order to derive the pre-exponential factors and determine 
the asymptotic behavior of P(n). We will now show how to obtain the results 
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for each of the regions. 



3.3.1 The non-critical regions (S± and So) 
3.3.1.1 £_ (h < -1) 

In this region (7 7^ 0, h < — 1) the generating function (13. 9p is nonzero for 
all q (see Fig. 13.2b ): this is the simplest case and can be treated using the 
(rigorously proven) Strong Szego Limit Theorem, see (1B.5j) . It gives 



P{n) = |det(S n )r~ 00 J B„(/ i , 7 )e- /3(h ' 7)n (3.40) 
with /3(h, 7) given by ([312]) or <KWi and 

{h, 7) = exp ( J2 ka k a- k j , (3.41) 



>fc=i 



where &k is defined in ( 1B.7[) as the fc-th Fourier component of the logarithm 
of a: 



°k = [ p- [logcr(g)]e" 
Jo 27T 

log [ 14 



27r dq e _ ifc?1 __ I , , cos g - /i + i 7 sing 
^ n V \/ ( cos 1 ~ h) z + 7 2 sin 2 g 



The sum in (13.411) is convergent only for 7 7^ and for h < —1. For h > 
— 1, the generating function (13.91) develops singularities which produce 1/k 
contributions to (I3.42p that make the sum in (13.411) divergent. Therefore, in 
the rest of the phase diagram these singularities have to be treated to absorb 
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the harmonic series contributions. Consequently, each region of the phase 
diagram will involve a different definition for the pre-exponential factor and 
the " regularization" procedure will sometimes generate an additional power- 
law contribution. The result is given by the Fisher-Hartwig conjecture that 
we must use in the remainder of the phase diagram. 

3.3.1.2 S (\h\ < 1) 

As can be noticed from Fig. 13.2b . in So (7 7^ 0, —1 < h < 1) the generating 
function a(q) vanishes and its phase has a discontinuity of tt at q = n. The 
asymptotic behavior of Toeplitz determinants with this type of singularities in 
the generating function is given by FH, which is actually proven for cases in 
which only one singularity is present. 

We decompose the generating function as in ( IB. 81) 



a(q) = r ( ? ) e ife-T)-od2 7 r- 7 r] (2 _ 2cos(g - tt)) 1/2 (3.43) 
and using (IB .91) we obtain 

P(n) = |det(S n )| n ~ oo J E (/i, 7 ) e -^ n . (3.44) 



The behavior is exponential as before with the decay rate f3(h, 7) from 
but the pre-exponential factor is different. According to (IB.10I) it is given by 



^,7) = ^ (3-45) 
t_(tt) 
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where, as in (IB. 6ft and (IB. 71) 

E[t\ = exp ( kf kr-k ) (3.46) 



oo 



,fc=l 



and 



% = ° k - 0(k). (3.47) 

Here 8(k) is the usual Heaviside step function. As we mentioned in the previous 
section, &k (I3.42[) has 1/k contributions from singularities of a(q) and the 
effect of the parametrization (13.431) is to cure (remove) these harmonic series 
divergences of the pref actor of the EFP in this regime. 

3.3.1.3 £+ (h > 1) 

In S + (7 ^ 0, h > 1), a(q) vanishes at q = and q = it and its phase presents 
two 7r jumps at those points (Fig. 13.2b ) . 

In this case the application of FH leads to some ambiguity, because there 
exist two representations of the kind ( 1B.8I) and one obtains two values for 
P(h, 7) using the two representations of the generating function: (3\ = f3 and 
P2 = P + iTTj with (3 from (13.121) or (I3.33[) . This ambiguity is resolved by the 
(yet unproven) generalized Fisher-Hartwig conjecture (gFH), which gives EFP 
as a sum of two terms so that both values of /?'s are used (see in the Appendix 
[R3lor [37]). 

The two leading inequivalent parametrizations (1B.15I) are: 
a(q) = r i( 9 ) e ^[(9-T)mod2^ ] (2 _ 2 cos ( ? _ 7r )) 1/2 

xe -l[ 9 mod2 7 r- 7 r] ^ _ 2cosq ) 1 / 2 (3.48) 
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T 



2( g ) e -i[(?— ) mod 2™] ( 2 - 2cos(g - tt)) 1/2 



X 



e i[gmod2^]( 2 _ 2cosg )l/2_ 



Application of flB.16j) gives the asymptotic behavior of the determinants as 



|det(S n )| n ~°° [E*(h t >y) + (-l) n E*(h n )] e~^ n (3.50) 



with 



^.7)^^L 0.52) 



and j3(h, 7), -E[r] defined in ( I3.12|I3.46I) with 



t k = a k - 9{k) - 1 9{—k). (3.53) 



Once again, as in the previous section, the effect of the parametrization is 
to remove the 1/k contributions to at ( 13.421) due to the singularities of the 
generating function. 

We conclude that the non-critical theory presents an exponential asymp- 
totic behavior of the EFP. In the region S + , however, the EFP in addition has 
even-odd oscillations 

P(n) n ~°°El(h,7) [1 + A + (/i, 7 ) cos(Trn)] e"^'^, (3.54) 

where the exponential decay factor is given by (13 .33j) . 
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The amplitude of the oscillations is 

T + (0)T_(7T) 



A + (h n ) 



r_(0)T + (7T) 

r(0) /t_(tt) 



t(tt) \r_(0) 



where we used ( IB. 1 1|) . the definition of r and (1B.13j) . We can deform the 



contour of integration as in Fig. I3.4f and calculate the integral in (13.551) to 
obtain 



A + (h, 7 ) = Vmn) = — . ( 3 - 56 ) 

+ Y - 1 + 7 

where K{h,^) was defined in (13. 2 1 p . 

Expression (13.541) for the EFP fits the numerical data remarkably well 
(see Fig. 13.51) and this fact strongly supports the generalized Fisher-Hartwig 
conjecture. 

One can understand these oscillations as a result of "superconducting" 
correlations of real fermions described by the Hamiltonian (12.61) . Fermions are 
created and destroyed in pairs of nearest neighbors. At large magnetic fields, 
the oscillations are due to the fact that the probability of having a depletion 
string of length 2k — 1 or 2k is very similar. Since the magnetic field in (12.61) 
is essentially a chemical potential for the fermions, the energy cost to destroy 
a pair of particles is Ah: at very big magnetic fields, the amplitude for a pair 

2 

destruction event is suppressed by a factor of ^, i.e. a probability of 
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This means that the probability of depletion behaves like: 



— J and 

UhY 2k 



P{2k) ~ l-j , (3.57) 



where the factor of two in the first expression is a simple combinatorial factor. 
The two probabilities in ( 13. 5 7ft can be combined in a single expression: 



P{n) = E [1 + A cos(Trn)] ( — J , (3.58) 



which is precisely (I3.54[) . with 



h \h 2 

We can check the correctness of this interpretation by taking the limit of 
(13T54D for h » 1,7. From (13321) and (l3~56l) it is easy to find 



^-00,7) = log^+oQ^) (3.60) 
A+(/i-oo,7) = 1 -X + (/?) (3 ' 61) 
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Figure 3.5: Results of the numerical calculation of the Toeplitz determinant are 
shown as points, P(n)e^ n vs. n at 7 = 1, h = 1.1. The value of (3 is obtained 
numerically from (13.351) . The solid line is the analytic result E(l + (— l) n A) 
with A = 0.2182... from (13.561) and E = 0.6659... obtained by fitting at large 
n. To make the plot more readable we show only every 11th point (for n = 
1, 12, 23, . . .) of the numerical calculation of the determinant. Note that the 
size of the points is not related to the estimated error in the numerics, which 
is actually smaller. 
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3.3.2 The critical lines (Q±) 
3.3.2.1 tt + (h = 1) 

For h = 1 the generating function a(q) vanishes at q = ir and its phase has 
7r jumps at q = 0, ir (see Fig. 13. 2H ). As in the previous section, the existence 
of two singular points gives rise to many terms of the form ( IB. 151) . However, 
in contrast to the E + region, the application of gFH as in (IB. 161) shows that 
all terms are suppressed by power law factors of n with respect to the leading 
one. 

The leading term is generated by the parametrization: 



a 



(q) = T 1 {q)e^ [(q - W ) mod2 — 1 (2 - 2cos(g - tt)) 1/2 e -3[<? mod2 ---] (3.62) 



and consists of an exponential decay with (3(1, 7) from (I3.35j) and a power law 
contribution with critical exponent A = 

|det(S n )| ~ Ej^n-Te e -W,7)n (3 63 ) 



with 

E\ (7) = E[t]G (f) G (?) , (3.64) 

W VV 2 1 /4 r V 4 (o) r _(7r) 

where G is the Barnes G-f unction defined in (IB.12j) and E[t] is defined as in 

(B36D with 

1 , ^ fc \ 1 1 



h = o k + ( - - (-l)M - 6(k) - - 9(-k), (3.65) 



with a k from (l3T42|) . 
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The next term (subleading at n — > oo) is obtained from the parametrization 
a(q) =T 2( ? ) e -i[(3- 7 r)mod27r-7r] ( 2 - 2 cos(g - tt)) 1/2 e 1 ^ 9 mod (3.66) 

and is given by 

Ef(-y)(-l) n nT^e-^ n (3.67) 

with 

E! (7) S E[r]G (I) G (I) £<°> , - . (3.68) 



M ' W 2 3 / 4 t! /4 (0)t + (it) 
Although the inclusion of the latter (subleading) term is somewhat beyond 
even gFH, we write the sum of these two terms as a conjecture for EFP at 
h = 1 

P(n) ~ Ej (7) nT fl+ (-l) n A 1 ( 7 )/n5j e ^ (1 ' 7)n . (3.69) 

As these results rely on our unproven conjecture, we present our numerical 
data for this case in Fig. 13.61 Indeed, we see that the form (13.691) is in good 
agreement with numerics and this supports our hypothesis. 

The amplitude of the oscillations is 

AM = 1_ G (l) G (l) r + (0)r_(7T) 
1 r(|) r(0) fr^Y 



V2 T(i) t(tt) Vt-(O) 

r (!) 1 (av I dz lo S^) 

exp Ulim*- — „ .„ , (3.70) 



r(i) 7 * V 7 271 z 2 - (1 + e 
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n 



Figure 3.6: Results of the numerical calculation of the Toeplitz determinant 
are shown as points, P{n)e^ n n l / l& vs. n at 7 = 1, h = 1. The value j3 = 
log 2 + 2G/ir with Catalan's constant G is obtained from (13.121) . The solid 
line is the analytic result E(l + (-l) n A/n^) with A = 0.2399... from fl3T2l 
and E = 0.8065... as obtained by fitting at large n. To make the plot more 
readable we show only every 11th point (for n = 1,12, 23, . . .) of the numerical 
results on the determinant. Note that the size of the points is not related to 
the estimated error in the numerics, which is actually smaller. 
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where we used ( IB.llj) and the identity 



G(z+l)=T(z)G(z). (3.71) 

To calculate the integral we deform the contour of integration as in Fig. I3.4fe 
and find 

r (?) 1 

Alfa) = — )# -f=- (3-72) 

iyv r(|) V27 1 ; 

We conclude that at h — 1 the EFP decays exponentially at n — > oo but 
with an additional power law pre- factor and a damped oscillatory component. 



Remark. It is curious to notice that the exponents 1/16 and 9/16 in (IX 

n 

and (13.671) remind us of the scaling dimensions of spins a x and a y . □ It looks 
as if the EFP operator (13.11) . among other things, has inserted square roots of 
spins transverse to the magnetic field at the ends of the string. 

3.3.2.2 (h = -1) 

For h — — 1 the generating function a(q) does not vanish but has a phase 
discontinuity of ir at q = ir. We parametrize a(q) as 

a(q) = r i( g ) e -ito— )mod2^^] (3 73 ) 

and apply FH to obtain 

P(n) ~ E 1 _ 1 fa) e-Pfa 1 '^ 71 (3.74) 



3 See Ref. [31] or (|2.43I2.45|) . where it was shown that the power laws for the a x and a v 
correlators are 1/4 and 9/4 respectively. 
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with 



5\ r m M 



where (3{— 1,7) and -Efr] are defined in (13.371) and (13.461) with 



(3.75) 



-iy 



4A- 



4A- 



6{-k) 



(3.76) 



and (Tfc from (13.421) . 

We can stretch the gFH the same way as in the previous section for h = +1 
by considering the second parametrization 



a(q) = r 2( g ) e ii[(<^)mod2^] 



(3.77) 



which gives 



P'{n) ~ El^i) n"n e 



A -/3(-l,7)« 



(3.78) 



with 

(i) a (I) |H 

Adding this subleading term to (I3.74p we obtain 



(3.79) 



P(ri) ~ #^(7) l + A_i(7)/n5 



-n/3(-l,7) 



(3.80) 



with 



A-i(7) 



GQ)G(J) r + (vr) 



G(f)G(|) r_(vr) 

r(|) r + (7r) 
r(l) r_(7r)- 



(3.81) 
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We propose (13.801) as an asymptotic form for EFP at h = 



-1. 



3.4 The line Fe- an exact calculation 



Before we conclude our analysis of the EFP with the study of the isotropic 
XY model, let us check our results (I3.40f3.12l) on the special linqj in the phase 
diagram defined by 

/i 2 + 7 2 = l. (3.82) 

It was shown in Ref. [M] that on this line the ground state is a product of 
single spin states and is given by 



iG>=niM=n[ c ° s (0 



cos(-l|U) + (-l) J sin(-)|U) 



(3.83) 



where | is an up-spin state at the lattice site j, etc. One can directly 
check that the state (13.831) is an eigenstate of (12. ip if the value of parameter 8 



is 



cos 2 6 



1-7 
1+7 



(3.84) 



and (13.821) is satisfied. It is also easy to show [31] that this state is, in fact, 
the ground state of (12. ip . 

The probability of formation of a ferromagnetic string in the state (13.831) 
is obviously 

P(n) = sin 2 " ( (3.85) 



2 2 \h\ V 1 + 7 



4 We are grateful to Fabian Essler who suggested us to check our results on this special 
line and pointed out the reference [33] to us. 
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which is an exact result on the line (13.821) . The value of /3(h, 7) which imme- 
diately follows from this exact result is 

(3(h = iv/1^,7) = - log Q =F = ~Hh,l), (3.86) 

where A(/i,7) was defined in (13.191) . 

This is, indeed, consistent with (13.231) since under the condition (13.821) the 
function A(/i, 7) vanishes. The integral (I3.20j) defining A(/i,7) vanishes for 
( I3.82p because the branching points (13 . 161) and ( 13. 17ft collapse to the same point 
and therefore the region of integration shrinks to just one point (I3.20p . In fact, 
the Toeplitz matrix ( 13.81) generated by ( 13.91) becomes triangular on the line 
( I3.82p with diagonal matrix element {S n )jj = sin 2 (6 l /2) and the determinant 
of S n is exactly (I3.85[) . 

From the definitions of /3(h, , ~f), we see that the decay factor consists of 
two terms, which now have a clear physical meaning. The A(h, 7) term is 
the factor we found above in (I3.86P and represents the contribution given by 
un-entangled spins. The remaining part accounts for the correlations between 
spins. Both A(/i,7) and the correlation functions given by (I2.22p and (12.231) 
vanish on the line (13.821) . 

Finally, it should be noted that there are actually two ground states for 
the theory on this special line: 

\G + ) = ]J\9,j) 
3 

\G-) = ]J\e,j) 
3 



n 

3 

n 



cosi- 



cosi- 



T,j) + (-l) J sin(-)||,j) 



T,J) 



■l) J sin(-)||,j) 



(3.87) 
(3.88) 
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These two states break the translational symmetry and are orthogonal in the 
thermodynamic limit. The reason for this degeneracy of the ground state is 
yet not well understood. 

If we were to consider a linear combination of these states, any local op- 
erator (involving only a finite number of lattice sites) would have vanishing 
cross-terms and very likely the same expectation value on either state (though 
one can design an operator which would violate the latter property). There- 
fore, for our present interest, evaluation of the EFP on any combination of 
these two ground states would yield the same result (13.851) and we don't need 
to concern ourselves with this degeneracy. 

The existence of these two ground states is instead of great interest in the 
study of quantum information applied to the XY model. We recently noticed 
this degeneracy and this helped us understand a previously known formula 
on the entropy of a block of neighboring spins [35]. This entropy is known to 
quantify the degree of entanglement of two subsystems of a system (in this 
case the block of neighboring spins and the rest of the chain). We remark 
in passing that many of the approaches to the calculation of the entropy for 
the XY model use Toeplitz determinant representations and the theory on 
Toeplitz determinant for deriving their asymptotic behavior, in the spirit of 
the present work. 
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3.5 The critical line (7 = 0) and the Gaus- 
sian behavior 



The case 7 = 0, corresponding to the Isotropic XY Model, has been studied in 
Ref. [24]. For 7 = the generating function (13.91) is reduced to the one found 
in [21]. 

For \h\ < 1, the generating function a(q) has a limited support between 
[— cos" 1 h, cos" 1 h]. To find the asymptotic behavior of the determinant of the 
Toeplitz matrix one can apply Widom's Theorem [38] and obtain [23] 

P(n) ~ 2A e ^'(" 1 )(l - ^-sn-i 2 • ( 3 - 89 ) 

We see that in this case, the EFP decays as a Gaussian with an additional 
power-law pre-factor. 

In a different context, the formula (13.891) appeared also in [5S] as a proba- 
bility of forming a gap in the spectrum of unitary random matrices. This is not 
unexpected, since the joint eigenvalue distribution of unitary random matrices 
is known to coincide with the distribution of free fermions in the ground state. 

For \h\ > 1, the theory is no longer critical and the ground state is com- 
pletely polarized in the z direction, giving a trivial EFP P(n) = for h > 1 
and P(n) = 1 for h < — 1. 
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3.6 Crossover between Gaussian and exponen- 
tial behavior: a Bosonization approach 

In order to understand qualitatively the crossover between the Gaussian asymp- 
totic behavior at 7 = and the exponential decay for 7 7^ 0, we employ a 
bosonization approach similar to the one used in [23], which will be sketched 
in Chapter In the limit 7 < 1 we consider the continuum limit of (12. 7p . 
bosonize the fermionic fields, and write the Euclidean action of the theory as 
S = J dx drC, where r = it is the imaginary time and the Lagrangian is 



where we also rescaled the time r — > VpT, with vp = 2\/l — h 2 , the Fermi 
velocity at 7 = 0. 

This is a Sine-Gordon theory for the "conjugate field" i?(x,t), which de- 
scribes the imaginary time dynamics of our 1-D system. In terms of d the 
density of fermions is given by p = A^d T '& + p , where p = k F /ii is the density 
of fermions in the ground state. 

In the field theory approach, the EFP (see Ref. [23]) in the limit n — »• 00 
would be given with exponential accuracy by the probability of an instanton 
P(n) ~ e~ So , where So is the action of the instanton. Here the instanton is 
the solution of the classical equations of motion of (I3.90p which corresponds to 
the formation of an emptiness of length n at the time r = 0. Unfortunately, 
the EFP instanton involves large deviations of the density of fermions from 
the equilibrium density po and is beyond the bosonization approach as the 




(3.90) 



63 



derivation of f 1 3 . 9 j) relies on the linearization of the fermionic spectrum near 
the Fermi points. 

Following [23], we are going to slightly generalize our problem by consid- 
ering the depletion formation probability instead of the EFP requiring 

PUo, o<x<n = Po + ^ 9 t#(x, r)\ T=0> 0<x<n = Po-p, (3.91) 

where p is some constant. The original EFP problem corresponds to p — po- 
Here, instead, we consider the probability of weak depletion, i.e. 

" ^ d T #(x, r) | r=0 _ 0<x<n = p«p . (3.92) 

We study the latter using an instanton approach to ( 13.901) and infer the (qual- 
itative) behavior of the original EFP from this weak limit. 

To simplify the problem further, we assume that the instanton configura- 
tion is completely confined to one of the wells of the Cosine potential in ( 13.901) 
and that the field d is small enough to allow for an expansion of the Cosine: 

S « ^ J dx dr [{d^f + 4 7 tf 2 ] . (3.93) 

In this formulation, the anisotropy parameter 7 1 / 2 plays the role of the mass 
of the bosonized theory. The probability we are looking for is given by the 
action Sq of the classical field configuration which satisfies the Euler-Lagrange 
equation (in this Klein-Gordon equation in two dimensions) with the 

boundary condition ( 13.911) 

p-(n) = e~ So . (3.94) 
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In the limit 7 = 0, the theory is massless and scale invariant. In [23] 
it was shown that, due to the scale invariance, the action of the instanton 
is quadratic in n. The instanton configuration in this case is essentially a 
droplet of depletion in space-time with dimensions proportional to n both in 
the space and time direction, in order to satisfy the boundary condition (I3.9ip . 
This result is consistent with the Gaussian asymptotic behavior prescribed by 
Widom's theorem (see Sec. 13.51) . 

In the massive case, a finite correlation length £ ~ <-y -1 / 2 is generated 
and we observe a crossover. For string lengths n smaller than the correlation 
length 7 -1 / 2 , the instanton action is not sensitive to the presence of the finite 
correlation length and is still quadratic in n (giving a Gaussian decay for EFP). 
In the asymptotic limit of string lengths greater than 7 _1//2 , the time dimension 
of a depletion droplet is of the order of £ (instead of n as in the massless limit): 
the action is linear in n and the probability has an exponential behaviorlfl 

In the next section we show how to solve the integral equation correspond- 
ing to the boundary problem (13. 91113. 931) and present its numerical solution 
and some analytical results. Figures 13.71 and 13.81 clearly show the crossover 
between a quadratic behavior of the stationary action for small n to a linear 

asymptotic one for n — > oo. 

5 This picture is very similar to the one for massless theory at finite temperature. In the 
latter the inverse temperature plays the role of the correlation length [53] fsee !3.7p . 
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n 



Figure 3.7: Plot of the value of the Figure 3.8: Plot of the derivative 
stationary action S vs. the string dS /dn with <S from (I3.103p . The 
length n. The action S is obtained plot corresponds to m = 0.01, p = 0.2 
from (13.1031) with f(y) given by the and clearly shows a crossover from the 
numerical solution of the singular in- quadratic to the linear behavior at 
tegral equation f)3.10ip . The graph n ~ 2/m = yl/7 = 200. 
depicts So(n) for m = 2^/7 = 0.01, 
p = 0.2. The crossover takes place 
around n ~ 2/m = 077 = 200 - 

3.7 Calculation of the stationary action in the 
bosonization approach 



In the previous section we have formulated the XY model near 7 = in 
terms of the bosonic field with Lagrangian (I3.93p . It was also pointed out 
that, instead of the EFP, we are interested in the Probability of Formation of 
Weakly Ferromagnetic Strings (PFWFS) and that we are going to calculate 
this probability in the saddle point approximation. In this section we are going 
to show how to solve the integral equation that solves this problem. 

This is a fairly technical section and can be easily skipped if the reader is 
not interested in the mathematical techniques required to tackle this problem. 

We consider a configuration of the field (instanton) which satisfies the 
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d T , d(x, r)| 



r=0,0<x<n 



(3.95) 



and that minimizes the action, i.e. that satisfies the Euler-Lagrange equations: 



The latter equation is the Klein-Gordon equation with the mass given by 
m 2 = 47 (see (13.931) ). The PFWFS will be found from the value of the action 
So corresponding to this instanton configuration ( 13.941) . In this appendix we 
calculate the stationary action needed in Sec. 13.61 

We now solve the differential equation (13.961) with non-trivial boundary 
condition (I3.95P by recasting it as the integral equation: 



where K (x, x'; r, t') is the modified Bessel function of 0-th order - the kernel 
of the differential operator (I3.96P in two dimensions. We impose the boundary 
condition (13.95}) by requiring that the "source" f(y) satisfies 



(d^ - m 2 ) $ = 0. 



(3.96) 




(3.97) 




y/np. 



(3.98) 
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This is the integral equation on f(y) we have to solve. 

Once the limit r — > is taken, the kernel in Eq. ( 13.981) is singular. We 
isolate the singularity by rewriting equation (I3.98P as: 



d i r f(v) i r 

1 lS0L d y + lim G (x,T ]y )f(y)dy = 2V^P (3.99) 

dx 7T Jo x - y t~>0 ir J 



with 



or, after integration over x, as 

i r f(v) i r 

/ ±S0L dy + _ G(x;y)f(y)dy = 2V^px 
K J Q x-y re J 



(3.100) 



(3.101) 



with 



G(x; y) = lim / G ( Xl , r; y) da*. (3.102) 

We have recasted Eq. (I3.98P in the standard form for a singular integral equa- 
tion (I3.10ip . Once we have the solution of this equation, we can calculate the 
action corresponding to this instanton as 



S = ^£ f(y)dy. (3.103) 



We solved the singular integral equation (13.1011) numerically and we com- 
puted the corresponding action (13.1031) . The results of this calculation are 
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presented as a plot of the action <S vs. n in Fig. 13. 7\ where we notice the 
crossover from a quadratic to a linear behavior (corresponding to a crossover 
from Gaussian to exponential behavior for the probability, (13.941) ) as we ex- 
pected. To confirm the nature of this crossover, in Fig. 13. 81 we plot dSo/dn and 
we see that it starts linearly and then saturates asymptotically as it should. 
In the limit n <C 1/m, we can expand the Bessel functions in the kernel 



G (x,r;y) = 



m 2 f T 2 1 , r, .o 21 i m ^ 1\ 

( ^7 + i h [('-») 2 + T*] + ln- + G--j + ... 

(3.104) 

where G is Catalan's constant. Then we solve the singular integral equation 
( 13.1011) to first order by first transforming it into a regular integral equation. 

In [39], Chap. 14, Sec. 114 it is explained that a singular integral equation 
like ( 13.10ip is equivalent to 

1 f n 

f{x) + - / N(x; y) f(y) dy = 2py/nx(n - x), (3.105) 

7T1 







where the new kernel is 



_ y/x{n-x) [ n G{y'-y) , 

N(x;y) = f -dy. (3.106) 

vri Jo ^ y '{n-y')(x - y') 

Using (13.1041) . we can explicitly calculate the integral defining N(x;y) in 
terms of elementary functions and after some algebra the integral equation 
(13.1051) results in a long, but essentially simple, regular integral equation. Its 
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n 



Figure 3.9: The solid line is the plot of the stationary action ( 13.1081) against 
n. This analytical solution is valid for n « 1/m and corresponds to m = 
0.01 and p = 0.2. The dotted line represents the value of the action (13.1031) 
with the source given by numerical solution of the singular integral equation 
(13.1011) . The dashed line corresponds to the zeroth-order, pure Gaussian, 
solution, i.e. (13.1081) with m = 0, which we include for comparison. We 
see that the inclusion of first order correction almost doubles the range in 
which the analytical solution is accurate. 
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solution is 



/(*) 



\ m 2 n 2 ( mn 3\1 r - 



n — x) 



/—_ m 2 n 2 ( n\ _ 1 / 
-VTrp — ^— (z- 2 J tan y 



(3.107) 



The corresponding stationary action (13.1031) is 



2 irp 



2^2 r 



2 2 
m n 



( ln _ + G _ 2 ) 



S = n 



8 



1 + 



16 



(3.108) 



The first term in (13.1081) corresponds to the Gaussian decay of PFWFS we 
expect in the limit of m = 0. In Fig. 13.91 we compare this analytical result for 
the action with the numerical result of Fig J3.71 In the plot, we include the pure 
Gaussian decay (the first term in (13.1081) ). which already gives a remarkable 
agreement for small n. The full solution (13.1081) extends this agreement further 
for larger n. For m = (7 = 0), (13.1081) reproduces the result calculated in 
[23] for h = 0. 

3.8 Emptiness Formation Probability at finite 
temperature 

Finally, we consider what happens at finite temperature (T > 0). The corre- 
lators and (127231) become 
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Gj k = (^I)t= r^^(l + cos^tanh^)e^- fc ), (3.109) 

J 

where the additional factor takes care of the thermal average. 
The EFP is expressed by 




(3.110) 



and in the spinless fermion formalism it becomes 

n 

P(n) = (H^I)t. (3.111) 

i=i 

We again use Wick's Theorem (or its thermal version, called Bloch-de 
Dominicis theorem [40J) to express it as a Pfaffian. The calculation proceeds 
the same way as for zero temperature and the EFP can be represented as 

P{n) = |det(T n )|, (3.112) 

where T n is the n x n Toeplitz matrix generated by the function 

t(q) = \ (l + eP'tanh , (3.113) 

where the "rotation angle" and the spectrum e q were defined in (12.91) and 
(12. lip respectively. 

The generating function t(q) is never-vanishing and has zero winding number. 
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Therefore, for T > we can apply the standard Szego Theorem to obtain 



P{n) n ~°°E{h,j,T)e 



-n/3(h,y,T) 



(3.114) 



where 



with 



/?(/*, 7, T) 



r% log I* Mi 



2?r 



dq 
2^ 



log 



^l + ^^tanh^ 



E(h,j,T) = exp ^ fc£ fc £- 



,fe=i 



, (3.115) 
(3.116) 



2?r 



dg 
— < 

2vr 



"i*9 l c 



1 / cos g - /i + i7 sin g e„ 

1 H tanh — — 

£„ 2T 



(3.117) 



and e q is given as in (12.111) by 



J (cos q — h) 2 + 7 2 sin 2 g. 



(3.118) 



As can be expected from simple thermodynamic considerations, at finite 
temperature the behavior is always purely exponential asymptotically. As 
it was shown in [23], at finite but very low temperatures one can observe a 
crossover from the zero temperature behavior at short string lengths n to the 
exponential behavior (13.1141) in the limit of very large n. This crossover occurs 
at a length scale of the order of the inverse temperature. 
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EFP for the Anisotropic XY model 


Region 


7, h 


P(n) 


Eq. 


Section 


Theorem 


£_ 


h < -1 




|3.40| 


13.3.1.1 


Szego 




-1 < h < 1 


Ee~ n(j 


13.401 


13.3.1.21 


FH 


£+ 


/i > 1 


E[l + (-l)M] e - ^ 


13.541 


13.3.1.31 


gFH 


Is 


7 2 + /i 2 = 1 




13.851 




Exact 




/* = 1 


Err 1 ' 1 * 


1 + (-1)M/Vn] e - ^ 


3.69 


3.3.2.1 


gFH 


n 


h = -l 


En- 1 / 16 




13.801 


13.3.2.2 


gFH 


^0 


7 = 0, < 1 


E n -V4 e -^a 


13.891 


E5 


Widom 



Table 3.1: Asymptotic behavior of the EFP in different regimes. The exponential decay rate (3 is given by Eq. (13. 12ft 
for all regimes. The regions in boldface are the critical ones. The coefficients E, A are functions of h and 7 whose 
explicit expressions are provided in the text. Relevant theorems on Toeplitz determinants are collected in the IB1 



3.9 Discussion and conclusions 



The asymptotic behavior of the Emptiness Formation Probability P(n) as 
n — > oo for the Anisotropic XY model in a transverse magnetic field as a 
function of the anisotropy 7 and the magnetic field h has been studied. We 
have summarized our results in Table 13.81 These asymptotic behaviors were 
first presented in [30] . 

Our main motivation has been to study the relation between the criticality 
of the theory and the asymptotics of the EFP. Let us now consider the results 
on the critical lines (fi an d fi± in Fig. 12. ip . The Gaussian behavior on O 
(7 = 0, \h\ < 1) is in accord with the qualitative argument of Ref. [23] using 
a field theory approach. In So (7 7^ 0, \h\ < 1) the asymptotic decay is 
exponential. We proposed a physical interpretation of the crossover between 
the two asymptotes using a bosonization analysis of the region of small 7: we 
suggest that there is an intermediate regime of Gaussian decay for the string 
lengths smaller than 1 / ^Ff which crosses over to the exponential behavior for 
longer strings. 

On the critical lines Q±, the decay of the EFP is exponential instead of 
Gaussian, and apparently contradicts the qualitative picture of Ref. [23] . The 
reason for this disagreement is that although at h = ±1 the model can be 
rewritten in terms of massless quasiparticles x defined in (12. 8p . we are still 
interested in the EFP for the "original" Jordan- Wigner fermions ip. In terms 
of x this correlator has a complicated (nonlocal) expression very much different 
from the simple one (13. 2p . From the technical point of view, the difference is 
that in the qualitative argument in favor of a Gaussian decay of EFP for 
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critical systems there is an implicit assumption that the density of fermions 
(or magnetization) is related in a local way to the field responsible for the 
critical degrees of freedom (free boson field <p). This assumption is not valid 
on the lines h — ±1. The theory is critical on those lines and can be described 
by some free field <fi. However, the relation between the magnetization and 
this field is highly nonlocal and one can not apply the simple argument of [23] 
to the XY model at h — ±1. 

Although EFP at the critical magnetic field does not show a Gaussian 
behavior, there is an important difference between the asymptotic behavior 
of EFP on and off critical lines. Namely, a power-law pre-factor n~ x appears 
on all critical lines. For the XY model it is universal (i.e. A is constant on 
a given critical line) and takes values A = 1/4 for 7 = [21] and A = 1/16 
on the lines h — ±1. It would be interesting to understand which operators 
determine these particular "scaling dimensions" of the EFP (see the remark 
at the end of Section 13.3.2.11) . 

At h > 1 the use of gFH predicts even-odd oscillations of P(n). We 
compared the predicted oscillations to numerical calculations of Toeplitz de- 
terminants and found a very good agreement (see Figs. 13. 5113. 61) . We proposed 
a physical interpretation of the oscillations as coming from pair correlations of 
spins which can be clearly seen as superconducting correlations in the fermionic 
representation (12.61) . 

In some parts of the phase diagram (£+, Q±) we used the so-called gen- 
eralized Fisher-Hartwig conjecture [37] which is not yet proven. However, our 
numeric calculations support the analytical results (see Figures [331 and [376]) . 
We note that to the best of our knowledge this is the first physically motivated 



76 



example where the original Fisher-Hartwig conjecture fails and its extended 
version is necessary. |^| We also suggest that the gFH could be used to find the 
subleading corrections to the asymptotic behavior, as we did for h — ±1 in 



be interesting to confirm it analytically. 

In conclusion, we notice that it was straightforward to generalize our results 
for nonzero temperature. The only modification is that at T 7^ the thermal 
correlation functions must be used instead of (j2.22|2.23p . Then, the generating 




function ( 13.91) is non-singular everywhere and we have an exponential decay of 
P(n) in the whole phase diagram according to the standard Szego Theorem 
and standard statistical mechanics arguments. We presented results for T^O 
in Section I3~8l 



6 We note that recently the theory on Tocplitz determinants has been used and extended 
with new results in order to calculate yet one more important physical quantity. We refer 
the interested reader to [41], [42] and [43], where the entanglement for the XY Spin chain 
and for Random matrix models have been calculated. 




. This novel hypothesis is supported by our numerics and it would 
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Chapter 4 

The Hydrodynamic Approach 



Very often it is possible to treat a quantum many body system as a continuum. 
A description of this kind views the system as a fluid, where the motion of the 
individual particles is not important, but one is interested in their collective 
behavior. This idea has a long history (see [T71 HU H3 [19] for instance). 

We consider a one-dimensional system at zero temperature, in the ther- 
modynamic limit. We will take a "semi-classical" limit, bringing Planck's 
constant h to zero and the number of particles N to infinity in such a way 
that their product stays constant. This limit allows us to describe the system 
only in terms of its local density and velocity. 

These two fields are clearly not sufficient to describe all states of the whole 
system, but they will provide an accurate approximation of it for some sec- 
tors of the theory. For instance we cannot consider configurations where the 
particles occupy disconnected regions of the phase-space (see Figure l4~2i) . 

The hydrodynamic approach can be applied to exactly integrable systems, 
where the description in terms of just two fields is essentially exact, once we 
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limit our attention to the 'allowed' sectors of the theory. 

The traditional hydrodynamic approach developed in the 1960's |Sj ad- 
dresses finite temperature systems and studies the time evolution of con- 
served quantities in many particle models. We concentrate instead on zero- 
temperature dynamics. 

The conservation laws guarantee that a perturbation or a fluctuation in the 
density of a conserved quantity of the system will not disappear on a length 
scale of the order of the interparticle distance, but will diffuse on a macro- 
scopic length scale. The diffusion equation is the main dynamical equation 
describing the evolution of these quantities. This description relies heavily on 
non-equilibrium processes taking place at a microscopical level to average into 
a diffusive motion. Once the dynamical equation is established, it is relatively 
simple to calculate some asymptotics of correlation functions for the conserved 
quantities using statistical mechanics methods. 

It is clear from this brief explanation that the hydrodynamic description 
is valid in the long wavelength approximation, with wavelengths longer than 
any other length scale of the system, so that diffusion averages have meaning 
and one can use the continuum description of the system. 

We are taking a zero temperature limit, where quantum fluctuations take 
the role of thermal fluctuations. 

Hydrodynamic equations are non-linear differential equations and their so- 
lution and treatment is not easy. If one linearizes them, the resulting the- 
ory is what is known in studying one-dimensional models as "Bosonization". 
Bosonization is a renowned tool, with a long history, to calculate correla- 
tors for the strongly interacting one dimensional systems. However, there are 
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phenomena which cannot be captured by a linearized theory. For instance 
the calculation of the Emptiness Formation Probability is beyond the linear 
approximation. In the next chapter we show how one can calculate the EFP 
to leading order for a number of systems using the hydrodynamic approach. 

Before we can analyze any system in particular, in section |4~T1 we introduce 
some definitions and show the basic structure of the quantities of interest. 
Then, in Section 14.21 we start explaining our approach on the example of the 
simplest possible system: Free Fermions (FF). This will allow us to warm up 
slowly to Section 14.31 where we generalize the approach to other systems with 
a derivation based on a Lagrangian Formalism. In section 14.41 we show how a 
linearized hydrodynamics reduces to the standard bosonization. 

Finally, in section 1431 we show how the Bethe Ansatz technique can be used 
to derive the hydrodynamic Hamiltonian of an integrable system. The Bethe 
Ansatz allows one to calculate exactly, although implicitly, the wavefunctions 
of an integrable system. Even more importantly, it provides fundamental 
information on the thermodynamics of the theory and an easier access to the 
conserved quantities of the model. For a short summary of the Bethe Ansatz, 
we refer the reader to Appendix O 

In Appendix [D] we show that if we neglect gradient corrections in a Galilean 
invariant system, the resulting hydrodynamic theory possesses a double infi- 
nite series of conserved quantities. It is still unclear whether these conserved 
quantities are enough to make the theory integrable in the sense of Liouville 
(i.e. the phase-space can be exactly factorized into a series of tori described 
by the action-angle variables). The answer is probably negative. But what is 
even more puzzling is the underlying symmetry guaranteeing this abundance 
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of conserved quantities and its physical meaning. We suggest that origin of 
this quantities should be connected to that of free fermions system, but we are 
unable to prove this conjecture at this moment. 

Zero temperature hydrodynamics was first developed by Landau in |17j . 
The collective description approach was successfully used in the context of 
lower-dimensional string theories and two-dimensional quantum gravity |46j . 
Our treatment is based on the work of A.G. Abanov and collaborators and it 
was described in [T9] . 



4.1 Some preliminaries 

Let us start by introducing some definitions. 

We want to develop a hydrodynamic description of a system in terms of 
its local density 

N 

p{x) = s ^^m5{x — Xj) (4.1) 



and current 



N 



-'sE (4 - 2) 

j = i v j 



where the sums are performed over the positions Xj of every particle in the 
system and 

{A,B} = A-B + B-A (4.3) 
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is the anti- commutator. The velocity is then defined as 




) 



(4.4) 



From this point onward we will always intend the velocity and current to be 
properly symmetrized, but we will not write it explicitly for the sake of brevity 
(e.g. (14.41) will be written as v = j/p). 

In a Hamiltonian formalism the dynamics of the system is encoded in the 
Hamiltonian (expressed as a function of the hydrodynamic parameters p and 
v ) and in the commutation relation between p and v. 

A "natural" assumption for the form of a hydrodynamic Hamiltonian is 



where e(p) is the internal energy per particle as a function of the local density 
of particles. The integration measure is equivalent to a sum over all particles, 
while the first term is the traditional macroscopic kinetic term. The second 
term in the Hamiltonian is some sort of potential term that includes all the 
interaction, but also the effect of the quantum statistics. We are going to see 
that this form is indeed the correct one in the next sections. 

A very important point is how we take the thermodynamic limit. We are 
going to see later that the internal energy e(p) is often an essentially quantum 
quantity, which vanishes in the classical limit (h — > 0). In such case, we 
are interested in taking the number of particles to infinity in a way that the 
internal energy stays finite as p — > oo and h — > 0. We are going to show how 




(4.5) 
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to derive the Hamiltonian for different systems in the next sections. 

Assuming standard commutation relation between position and momen- 
tum, 

[xi,pj] = XiPj - pjXi = ihSij, (4.6) 



from their microscopical definitions (14.1114.21) we can calculate 



\p(z),j(y)] = ihp(y)d,J(x - y) (4.7) 

or, using (14.41) . 

[p(x),v(y)] = -m'(x-y), (4.8) 

where 5'(x) denotes the derivative of the Dirac's delta-function. Eq. (14.81) 
is the canonical commutation relation for hydrodynamic parameters (see, for 
instance, [47]) and, in connection with the Hamiltonian of the system (14. 5p . 
will specify the time evolution of p and v: 



-d x (pv), (4.9) 



v 2 



- d 4j + (p £ )p ) • ( 4 - 10 ) 
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4.2 The simplest example: Free Fermions 

Let us consider a one dimensional system of Free Fermions (FF), without 
internal degrees of freedom. The local density of particles is given by 



I (4-11) 



Let us assume that the system moves with some velocity v(x), which itself 
is a slow function of x. This means that the left and right Fermi points are 
also functions of x so that: 

, x f kRix) dk k R ( X ) ~ k L (x) 

p[x) = m — r = m . 4.12 

The momentum density of the system is then given by 

where we recognized that the momentum of the system is in fact its current 

j(x) = p(x) v(x). (4.14) 

From this equation and using (14.12114.131) . we express the local velocity of the 
system as 

= k R (x) + k L (x) 
v ' 2m K J 

Equations (I4.12f4.15p give us the hydrodynamic parameters of the free fermions 
fluid p and v in terms of the left and right Fermi points ki,,R- 



84 




Figure 4.1: Depiction of the typical free fermions phase space configuration 
we can describe in our hydrodynamic description. At zero temperature, the 
particles are confined within the right and left Fermi Points (dark area) and 
these move in space (and time, not depicted here). 
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We invert (14. 12114. 151) and get the right and left Fermi momenta as functions 
of the density and velocity at each point in space and time: 



fin 

kn,L{x, t) =m v(x, t)± — p(x, t). (4.16) 



This is one of the most important limitations of our approach. For the Fermi 
points to be functions of the coordinates, they have to be single valued. This 
means that we can only describe star-like configurations (a.k.a. quadratic 
profiles) that in phase-space look like in Figure 14.11 We cannot address states 



similar to the one in Figure I4.2[ where disconnected regions of phase-space are 
populated or where the world-line of ku } L goes back above or below itself. 
The Hamiltonian of a free fermions system is 

fk R 7 2 JjL J.3 _ 1.3 

H= / — = (4.17) 
J k 2m 2irn Ylixnm 

and, by using (I4.16I) . can be written in terms of the density and velocity of the 
system as 

*(») = *= *^ + ^(«). 

We are now in a position to discuss what kind of thermodynamic limit we are 
interested in. As we increase the number of particles, we let the density grow 
as well. As we approach the semi-classical limit with h — > 0, we let the density 
grow so that the product of the two quantity remains constant: 

Hp — > const (4.19) 
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Figure 4.2: Phase space depiction of a free fermions system. The hydrody- 
namic description can be applied only in star-like configurations like the one 
between points b and c. Between a and b and between c and d the world 
line of the right Fermi Points kn(x) is not following a quadratic profile and 
the description of the system in terms of just two fields (density and velocity) 
fails. 
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This means that the energy per particle remains finite and that the kinetic 
and potential terms (respectively the first and the second term in ( 14.18P ) have 
comparable magnitude. To make this point more apparent, let us rescale the 
velocity as 

v -> — v (4.20) 
n 

so that the Hamiltonian becomes 



m z \ 2 bm z 



From this point on, we will set constants H = m = 1, this means that all 
quantities are going to be expressed in units of h and m: 



p(x)v 2 (x) 



7T 



2 



H{x) = v ^ + '-p\x). (4.22) 



Note that the second term on the right hand side gives the internal energy 
of the system as e(p) = ^p 2 and that this contribution comes only from the 
Pauli principle, since we are describing free fermions. For this reason it is not 
correct to identify the internal energy with a pure potential term. 

We showed in the previous section that the microscopic description of the 
system fl4.lll4.2ll4.4p imposes the following commutation relation between den- 
sity and velocity: 

\p(x),v(v)] = -ti'(x-y). (4.23) 
Using this relation, the Hamilton equations give us the dynamics of the system: 

pt(x) = [H(x),p(y)] = -d x (pv), 
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v t (x) = [H(x),v(y)] = - 




(4.24) 



We recognize the first equation as the continuity equation that expresses par- 
ticle conservation and relates density and velocity of the system as dependent 
quantities 

d t p + d x (pv) = 0. (4.25) 
The second dynamical equation is: 

d t v + vd x v = —n 2 pd x p (4.26) 



and is known as the "Euler equation" in the classical theory of fluids. 

We have shown in the simple case of free fermions how to construct a 
collective description of the system in terms of local density and velocity, 
which represents its hydrodynamic description. 

It is worth noticing here that we derived (14.25114.261) essentially semi-classically. 
For the free fermion case (and only for this case), it can be shown [46] that a 
quantum treatment of the theory gives exactly the same results, i.e. one only 
needs to consider commutation relations instead of Poisson bracket, all the 
functions are promoted to be operators (we put "hats" on top of the density, 
velocity and Hamiltonian) and all equations have to be interpreted as operator 
equations. 

In the next section we generalize this construction to some interacting 
systems. 
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4.3 Lagrangian formulation of Hydrodynam- 
ics 



As we showed in section I3.6[ the leading behavior of the EFP can be easily 
calculated in the semiclassical approximation. To this end it is useful now to 
turn to a Lagrangian formalism to calculate the hydrodynamic action and the 
partition function of the quantum theory [19] . 
We consider the partition function 



of some field or collection of fields u(x, t) and where u = dtu. 

Eventually, in the next chapter, we are interested in calculating the EFP 
as a rare fluctuation in the equilibrium configuration, an "instanton" , and this 
configuration will take place in imaginary time. To do this, we will need to 
formulate the Euclidean theory with the Euclidean action 



J J-1/2T 

with periodic boundary conditions in the imaginary time r = it defined by the 
temperature of system T. 




(4.27) 



where S[u] is the Action 




(4.28) 




(4.29) 
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Let us first consider a system with Galilean invariancaj. The requirement 
of Galilean invariance restricts considerably the form of the lagrangian to 



2 

C(p,v) = ^ + pe(p) + ... (4.30) 



where the first term is the kinetic energy of the system and is the only one 
that can depend on the velocity and the second term is the internal energy of 
the fluid, to be determined by the fluid's equation of state. e(p) is the internal 
energy per particle at given density p. Other terms can be included in the 
lagrangian and were here denoted by the dots: these terms include spatial 
derivatives of density and velocity and can be neglected if density and velocity 
gradients are sufficiently small. 

Particle conservation requires the continuity equation 

d t p + d x j = (4.31) 

to be satisfied, where j = pv. We can interpret this equation as a constraint 
on the fields, relating p and v to each other. We can solve this constraint by 
introducing the "particle displacement field" u(x,t) such that 

p = p + d x u, 

J = -d t u, (4.32) 
where po is the equilibrium value of the density at infinity, so that the displace- 



1 Lattice systems are more complicated and have not been treated in this thesis for lack 
of time, but their analysis will be addressed shortly. 
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ment field can satisfy standard boundary conditions and vanish at infinity. 
(One could define the displacement field in terms of the microscopic theory 
like in O4.104.2p as 

TV 

u(x) = 9(x — Xj) — pox (4.33) 
i=i 

where the sum is over the Xj position of all particles in the system and 6{x) is 
the Heaviside Step-function.) 

One can then write the Lagrangian in terms of the displacement field: 



u 2 

C(u, u) = — ■ + (p + u x )e(p + u x ) (4.34) 

2(p + Ux) 



and look for the field configuration minimizing the action S[u}. We can write 
the Euler-Lagrange equation for this Lagrangian in terms of the physical field: 



d t v + vd x v = -d x d p [pe(p)] , (4.35) 

where p and v were defined in terms of u(x,t) in 14.321 This is the Euler 
equation for a one dimensional fluid jH| and reduces to (14.261) for the free 
fermions case (e(p) = ^p 2 )- 

We see that, once the internal energy e(p) in (14.301) is known from the 
equation of state, one has everything to calculate the Lagrangian of a one- 
dimensional Galilean invariant system. In section 14751 we are going to show how 
to compute the internal energy e(p) as a function of the density for integrable 
models, using the Bethe Ansatz solution. Alternatively, for non-integrable 
systems one can compute e(p) from numerics or from other phenomenological 
methods. 
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The problem is not equally straightforward for systems without Galilean 
invariance, like for lattice models. Nonetheless, a hydrodynamic description 
can be developed in such cases as well. In chapter [6] we will consider systems 
with more than one type of particles, namely fermions with spin, and their 
hydrodynamics will not be as simple as the one for Galilean invariant models. 

4.4 Bosonization as a linearized hydrodynam- 
ics 

If the deviations from the equilibrium state are small, we can expand the 
Lagrangian (14.341) around p = p and j = as 

C(u,u)~±(u 2 + v 2 M u 2 x ), (4.36) 

where we defined the sound velocity at equilibrium as 

v 2 s0 = pd 2 p (pe(p))\ p=po . (4.37) 

By scaling the time as v s ot — > t we can write the linearized Action as 

S~— [ & 2 x\(d,u)\ (4.38) 

where d 2 x = dxdt and p, — x,t. 

The linearized hydrodynamics is described by the familiar quadratic action 
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of bosonization, giving the Laplace equation 



Au = 



(4.39) 



as the equation of motion. 

4.5 Hydrodynamics from the Bethe Ansatz 

The Bethe Ansatz is a powerful tool to analyze integrable systems. For the 
reader unfamiliar with this technique, we review in Appendix O the basic 
formulas that we need to develop our hydrodynamic approach. In fact, once 
one is familiar with the machinery, its use in the hydrodynamic formalism is 
quite straightforward. 

The Bethe Ansatz solution is based on the distribution of the quasi-momenta 
r(q). This distribution is found as the solution of an integral equation known 
as the Bethe Equation ( 1C.23I) : 



where the Kernel K(q — q') encodes the interactions of the system. 

For the ground state of the system, the limits of integration in the integral 
equation are chosen to be symmetric, i.e. from —k to k. But we are interested 
in a different sector of the theory, a sector characterized by a total finite 
momentum of the syste 

2 One can think as the ground state of the theory as coming from a grand-canonical 
ensemble approach. One then constructs the partition function by adding to the Hamiltonian 
a vector potential coupled to the total momentum of the theory, in addition to the standard 




(4.40) 
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We are guided by our experience with the free fermions model in section 

EE 



and we assume asymmetric limits of integration for the Bethe equation: 
r(q) + j"" K{q - q')r(q')dq' = -k (4.41) 



The number of particles and momentum per unit length of the system are 
given by (1C.25|C.26|) . If we assume them to be space dependent, we can iden- 



tify them as the density and current we need for our hydrodynamic description: 

rk R (x) 

p(x) = / r(q)dq, (4.42) 

Jk L (x) 

j(x) = / qr(q)dq. (4.43) 

Jk L (x) 

The energy of the system is given by 

rk R (x) 2 

H(x) = / q -r{q)dq. (4.44) 

Jkr.(x) Z 



This is an implicit expression for the Hamiltonian, since it depends on the 
momentum density r(q), which in turn is a function of kn and k^. 

We could invert (I4.42f4.43p and express kn and fci as a function of p and 
j (as we did for the free fermions) and then use this relation to express the 
Hamiltonian as a function of the hydrodynamic parameters density and cur- 
rent. 

We can actually do better. For a system like this we can use Galilean 



chemical potential coupled to the number of particles. 

3 One can derive the free fermions case from the general formalism we are developing in 
this section using the Bethe Ansatz. One has just to keep in mind that for free fermions, 
the kernel in (I4.40[) is equal to zero, therefore the momentum distribution is r = l/2n. 
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invariance to boost the reference frame. In practice, we perform a change of 
variable in the integrations by shifting the integration variable as 



where 



q = q-v (4.45) 



v = (4.46) 



so that equations ( 14. 42 [ 4. 43 1) become 

•k(x) 



p(x) = / T'(q')dq' (4.47) 

J-k(x) 

/+k(x) 
(q' + v) T'(q')dq' = p(x)v(x) (4.48) 
-k(x) 

where we use the fact the the function r'{q') = r(q + v) is even and where 



k = (4.49) 



More importantly, this change of variable writes the Hamiltonian as 



= ^ + P<P) (4.50) 



where 



1 f k a' 2 

e(p) = - / q -T>(q')dq'. (4.51) 
P J-k A 



This is precisely the expression we derived in generality in H4.5[) or (I4.30p . 
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but now we have a microscopic way to calculate the internal energy function 
e(p), from the Bethe Ansatz. Moreover, the formalism developed in this section 
will prove very powerful when applied to more complicated systems as we are 
going to see in Chapter EJ 

To conclude, we recall once more that from the commutation relation (14.81) 
and using (14.9114. 101) . we get the equations of motion (14.31114. 351) . 
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Chapter 5 



The EFP from Hydrodynamics 



We now turn back to the problem of calculating the correlator known as Empti- 
ness Formation Probability. In Chapter[3]we calculated the asymptotics of this 
correlation function in the XY model. This calculation was facilitated by the 
specific structure of the model which allowed us to express the EFP exactly 
as the determinant of a matrix. 

We argued in Chapter Q] for the importance of the EFP in the theory of 
integrable models. There, we also pointed out that the EFP can provide us 
with important insights in the general problem of calculating correlators that 
involve large deviations from the equilibrium configuration. We propose the 
hydrodynamic approach to address this problem. 

We introduced our hydrodynamic approach in the previous chapter; here 
we are going to show how this formalism helps us in calculating the leading 
asymptotic behavior of the EFP for some integrable models possessing a simple 
Galilean invariance. 

We are going to calculate the EFP as the probability of a rare fluctuation 
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of the theory, an "instanton" which depletes a region of particles. 

In section I57T1 we are going to explain how we are setting up the calculation 
and draw some general conclusion. In section 15.21 we are going to attempt 
the calculation using the linearized hydrodynamics (also know as Bosoniza- 
tion) and show that this approximation is only enough to produce the correct 
qualitative, but not quantitative result. In section 15.31 we are going to ma- 
nipulate the action of the theory to show that we can extract the leading 
behavior of the EFP from the asymptotic behavior of the instanton solution. 
In section 15.41 we are going to use this result to calculate the EFP for the free 
fermions hydrodynamics we constructed in section 14.21 Then, in section 15.51 
we will consider Calogero-Sutherland particles and calculate the EFP in the 
hydrodynamic approach. 

These results were presented for the first time in [19] . In the next chapters 
we are going to develop a hydrodynamic description for more complicated 
models and we will show for the first time how to calculate the EFP for these 
models. 

5.1 EFP as an instanton configuration 

We consider the Lagrangian formulation of the hydrodynamic theory intro- 
duced in section 14.31 We are going to work in Euclidean space, i.e. in imagi- 
nary time, because we are not studying the dynamical evolution of the system, 
but a property of the ground state. 
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We write the action of the Euclidean hydrodynamic theory as 

Z = J Vue~ SE[u \ (5.1) 

where Se[u] was defined in (I4.29p . 

Our approach is essentially the same as the one used in section 13.61 The 
asymptotic behavior of the EFP is defined as a rare fluctuation that depletes a 
region of length 2R from every particle. We interpret this configuration u(t, x) 
as the result of a collective motion where all particles move away so that at 
some time r = we have no particles in the spatial interval [-R, R]. 

This trajectory u(r, x) is a solution of the equations of motion satisfying 
the EFP boundary conditions: 

p(r = 0;-R<x< R) = 0, (5.2) 

and standard boundary conditions at infinity 

p -> p , x,r^oo, 

v — > 0, x, t— » oo, (5.3) 

Then with exponential accuracy, the EFP is calculated as 

P(R) ~ e~ Sopt , (5.4) 

where S opt is the value of the action (I4.29[) on the configuration u(x, t). 

In the introduction of Chapter [3] we already argued that we can estimate 
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the qualitative dependence on R of the stationary action S op t on very general 
grounds. The argument went as follows: in order to satisfy the EFP boundary 
conditions 15. 2[ the instanton solution needs to open a gap in its density that 
perturbs the configuration both in space and time. The spatial extent of this 
disturbance is clearly of the order of the gap to be opened, i.e. R. 

If we assume that the effective fluid we are describing is compressible and 
there are no other relevant length scales in the system, then the typical time 
scale of the disturbance has to be of the order of R/v s , where v s is the sound 
velocity at p = pq. Therefore, the (space-time) "area" of the disturbance is of 
the order of R 2 , the action S opt ~ R 2 and we conclude that the decay of the 
EFP (15.41) is going to be Gaussian: 



where a is some (non-universal) constant depending on the details of ( I4.29j) . 

In general, this behavior will persist until another length scale emerges to 
compete with R. At finite but sufficiently low temperatures, the temporal ex- 
tent of the instanton R/v s is smaller than the inverse temperature R/v s <C 1/T 
and the instanton does not feel the effect of the temperature, resulting in am 
intermediate Gaussian decay of EFP as in (15. 5p . However, as the temperature 
increases, or as we consider bigger R, the periodic boundary conditions in time 
become relevant and the 1/T scale defines the temporal size of the instanton, 
so that the space-time area of the disturbance scales as R and one obtains 



P(R) 



-aB? 



(5.5) 



P(R) ~ e 



-7-R 



(5.6) 
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Similarly, if some other term in the theory (14.291) drives the system away 
from criticality, the hydro dynamic approach will describe an incompressible 
fluid. The correlation length of the density fluctuations will have the same 
effect of a relevant length scale as just discussed and the EFP will decay 
exponentially as in (15.61) even at zero temperature if R is bigger than this 
correlation length. In both cases, as one increases R, one would observe a 
crossover from Gaussian to exponential behavior, with the crossover happening 
when Rr ^ v " or then R is comparable with the relevant correlation length. 

Once again, this is the same argument of [23] and we have shown the latter 
crossover phenomenon in ( 13. 6ft . 

Before we proceed to any actual calculation, for the sake of generality, 
we would like to consider a slightly different problem from the Emptiness 
Formation Probability (15. 2p . We are going to impose the Depletion Formation 
Probability (DFP) conditions on the instanton configuration: 

p(r = 0; -R < x < R) = p, (5.7) 

where p is some constant density. Setting p = we obtain the EFP problem 
while for p close to po one can use the bosonization technique to calculate 
P(R;p). 

In conclusion, we would like to find a solution to the equations of motion 
that satisfy the boundary conditions ( 15. 7115. 3ft and then use this solution to 
calculate the value of the stationary action S opt . In Euclidean space, due to 
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the different sign convention, the Euler equation is 



d T v + vd x v = d x d p [pe(p)], 



(5.8) 



and the continuity equation reads 



d T p + d x j = 0. 



(5.9) 



5.2 Linearized hydrodynamics or bosonization 

A simple qualitative analysis would tell us why we cannot calculate the EFP 
using a bosonization approach. In fact, one of the most important approxima- 
tions done in deriving the bosonization description of a system is to assume 
a linear spectrum. This is definitely a very reasonable assumption when one 
wants to look at low energy excitations close to the Fermi points, but the EFP 
clearly picks up contributions from every point in the spectrum, since it im- 
poses a strong deviation from the equilibrium distribution in the configuration 
space. As we pointed out in Chapter [U the inability of Bosonization to de- 
scribe the EFP is one of the main motivations for our study of this correlator, 
since it is a poster child of the advantage of using a hydrodynamic description 
over bosonization for certain calculations. 

However, we can try to calculate the DFP (15. 71) in the bosonization ap- 
proach if 



i.e., if we consider the probability of formation of a only small constant density 



Po - P 



< 1 



(5.10) 



Po 
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depletion along the long string —R < x < R. 

In this case the deviation of the density from po is small almost everywhere, 
as we will show below, and therefore we can assume that only particles close 
to the Fermi surface will be involved in the creation of the depletion and 
bosonization is applicable. 

As we showed in section 14 A\ we can derive bosonization by linearizing the 
hydrodynamic action (14.291) into ( 14. 38ft . The corresponding linearized equation 
of motion for the displacement field u(t,x) is the Laplace equation 

Au = (5.11) 

which has to be solved with the DFP boundary condition 

u(t = 0, x) — -(p - p)x, for - R< x < R. (5.12) 

In [19], this problem was solved by inspection and the correct instanton 
solution was found to be 

u(x, t) = -( Po -p)$l ^ - y/zi-R^j , (5.13) 

with the complex notation 

z = x + iv s0 T (5.14) 
and where v s o was defined in (14.371) as 

v 2 s0 = pdl(pe(p))\ p=p \l (5.15) 
1 It is important to remark that bosonization is universal and that it has the same form 
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One can check that indeed, at r = 0, —R < x < R the complex coordinate 
zq is real and square root in (15.131) is purely imaginary so that ( 15.12ft is satisfied. 
At space-time infinity Zq — > oo we have 



u(x,t) & — — — — (5.16) 



with 



a = a = - (p — p) R 2 (5-17) 



and where we used z = x — iv so r to denote complex conjugation. 

In terms of the original hydrodynamic parameters (14.321) . we obtain from 
( 15.161) that as zq — > oo 

a a 

Z Zq 

.v s0 ( a a\ 
v w -l — — - , (5.18) 
Po V^o zl) 

which obviously satisfy the boundary conditions at infinity (15.31) . 

By plugging this instanton solution (I5.13P into (14.381) . we can calculate the 
stationary action for the DFP problem as 



Sdfp = \— Kpo - p) Rf ■ (5.19) 
1 irpo 



Let us end this section by noticing that the instanton solution ( 15.131) is 
singular close to the ends of the string t = 0, x = ±R and therefore gra- 
dients of u, i.e. density and velocity, diverge. This is not consistent with 



for every theory. The specificity of each model is captured by the parameter v s q. 
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our approximation justifying the bosonization approach, i.e. that the solution 
does not deviate too much from its equilibrium value. However, in [19] the 
corrections coming from those areas were estimated to contribute only with 
terms of higher order in the small parameter (15.101) to the action (15.191) . 



5.3 EFP through the asymptotics of the solu- 



tion 



In [19], Abanov showed that the calculation of the value of the stationary 
action S op t ( 15.41) can be extracted from the asymptotics of the EFP solution 
of the hydrodynamic equations, using a Maupertui principle. 

The first step is to calculate the variation of the action (I4.29p with respect 
to the displacement field u 



5S P 



d 2 x \ — d T (vSu) — d x 



- d p (ep) Su 



-5u 



[d T v + vd x v - d x d p (ep)} J. 



(5.20) 



Note that we kept surface terms (full derivatives) in addition to the last 
term which produces the equation of motion (15.81) . 

Let us now calculate the derivative of the action Se with respect to the 
equilibrium background density po an d evaluate it on the EFP (or DFP) so- 
lution that saturates the equations of motion: 



pO <->opt 



d x < -d T (vd po u) - d x 



dp(ep) d po u 



(5.21) 
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We can rewrite the full derivative terms in (15.211) as a boundary term 



d P0 S O pt = j> vd po u dx + (d p {ep) - V — \ d po u dr 



(5.22) 



where the integral is taken over an infinitely large contour around the x — r 
plane. 

Equation (I5.22j) gives us the value of the (derivative of the) stationary 
action from boundary terms lying at infinitely large x and r. If we assume the 
asymptotic behavior for the instanton solution to be 



u(x, t) 



a a 



Po \ z z 

like in (15. 16115. 18p . after some algebra we conclude 



a a 

* « -i^h-l), (5.23) 



d P0 S opt = 27r^(a + «). (5.24) 
Po 



The problem of calculating the leading behavior of the EFP (DFP) is then 
reduced to the evaluation of the asymptotic behavior of the instanton solution. 
Once one knows a in (]5.23p . by integrating equation (I5.24p one finds the desired 
behavior of the EFP (DFP) problem. 

We can check this result with what we found in the previous section from 
bosonization. Substituting into (15.241) the value of a obtained in (I5.17P we 
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obtain 



d P0 S opt = ix— (po - p) R 2 (5.25) 
Po 



which is equivalent to ( 15.191) up to higher order terms in the small parameter 

(Em}. 

5.4 EFP for Free Fermions 

In section 14.21 we analyzed in details the hydrodynamic description of free 
fermions. The hydrodynamic equations for this system are 

d T p + d x (pv) = 0, 
d T v + vd x v = ir 2 pd x p. (5.26) 

with the sound velocity 



v s = v F = k F = pd 2 p (pe) = np (5.27) 

where we noticed that, since in our notation m — 1, the sound velocity and 
Fermi velocity are the same for free fermions. 
By introducing a complex field 

iv = np + iv, (5.28) 

we can rewrite both dynamical equations (I5.26P as the single complex Hopf 
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tiorfl 



equati 

d T w — iwd x w = 0. (5.29) 
The general solution of this equation can be written implicitly as 

w = F{x + iwr) = F(z), (5.30) 

where F(z) is an arbitrary analytic function of the complex variable z defined 



as 



z = x + iwr. (5.31) 

We can determine the function F(z) through the boundary conditions. In 
[19] it was shown that by defining 

F(*) = np + n(p - p) (5.32) 

equation (I5.30P satisfies the DFP problem and we can get the EFP solution 
by setting p — 0. 

One obtains the density and velocity by taking the real and imaginary part 
of w. To write down this instanton solution explicitly is quite complicated, 
since one has to solve the equation 



W -Txp = 7r(po -p)—j====. (5.33) 
sjz 1 — R 



Abanov in [19] plotted a numerical solution that we show in figures 15.11 and 



2 In real time formalism, instead of w and w, one introduces right and left Fermi momenta 
fcfi,L = irp±v which satisfy the Euler-Hopf equations dtk + kd x k = 0, reflecting the absence 
of interactions between fermions. 
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x_J 



Figure 5.1: The density profile p(x,r) is shown for the EFP instanton as 
implicitly given by f)5.33p . The density diverges at points (x, r) = (±i2, 0). 
The shape of the "Emptiness" is shown in figure 15.21 [From [19] ] 

However, as we argued in the previous section, to calculate the leading 
behavior of the EFP (DFP), we do not need to know the full solution of 
( 15.331) . but we just need the asymptotic behavior of the instanton solution 
using (I5.24p . 

In the limit x, r — > oo we have w — > np and z = x + iwr — > x + i7rpo r = z o- 
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Figure 5.2: The region of the x — r plane in which p(r,x) = for the EFP 
instanton solution for free fermions (I5.33l is shown. The boundary of the region 
can be found to be given by an astroid x 2 ! 3 + (npor) 2 ^ 3 = R 2 ^ 3 . [From [19J] 
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Therefore, w from (I5.33P behaves asymptotically as 



R 2 

w - np « 7r(p - p)—j (5.34) 

ZZn 



and taking, e.g., its real part we have 



J ^0 



Comparing (15.351) with (15.231) we identify 



a = - A {p Q -p)R 2 . (5.36) 



Substituting this into (15.241) we find 

d P0 S opt = 7T 2 (p - p)R 2 (5.37) 

and, after integrating in po, 

S op t = \ [Hpo - P)R? ■ (5.38) 
We can conclude that DFP and EFP are respectively 

P DFP (R) ~ exp|-^[7r(p -p)i?] 2 |, (5.39) 

Pefp(R) ~ exp|-^(7rp i?) 2 }. (5.40) 

We can compare this result with the known results on the EFP for free 
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fermions [36] that we recapitulated in Appendix |A] and confirm that (15.401) 
gives the correct first (Gaussian) term. 

5.5 EFP for the Calogero-Sutherland model 

Our final example for this chapter is the Calogero-Sutherland model, an inte- 
grable model of one-dimensional particles interacting with an inverse square 
potential. The Hamiltonian of the model is 



H 



j=l 3 l<j<k<N K 3 K > 

N 

- E 




A 



k=i,k& Xi Xk 



and we are going to consider it in the thermodynamic limit iV — > oo whil< 



keeping the density constant as we discussed in the previous chapter (J4T 



This model is known to be integrable [3H1 EH] and we briefly analyzed it 
in the appendix dedicated to the Bethe Ansatz flC.8jl . We can write explicitly 
the ground state wavefunction of (15.411) as 

^GS = H(x J -x k )\ (5.42) 

j<k 

As discussed in Appendix [Cj this formula shows that Calogero-Sutherland 



3 Technically, in order to go into the thermodynamic limit, we need either to add an 
external harmonic potential to confine the particles or to consider the model as defined on a 
closed ring; otherwise the repulsive nature of the interaction would spread the particles and 
we would not be able to keep the density fixed. We can neglect these details here, since they 
do not change anything relevant for us and, most of all, they both preserve the integrability 
of the model. 
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particles have intermediate statistics interpolating between non-interacting 
bosons (A = 0) and non-interacting fermions (A = 1). 

As a connection to the theory of Random matrices we mentioned in II .41 we 
should point out that that the probability distribution of particle coordinates 

i^Gsi 2 =ni x i-^i 2A ( 5 - 43 ) 

j<k 

at some particular values of the coupling constant A = 1/2, 1, 2 coincides with 
the joint probability of eigenvalues for the orthogonal, unitary, and symplectic 
random matrix ensembles respectively (see Eq. (11. 7p ). 

To calculate the leading behavior of the EFP for the Calogero-Sutherland 
model we need to determine the equation of motion of its hydrodynamic de- 
scription (15. 8p . which means that we need to know the internal energy function 
e(p). This can easily be found in [19] or from a direct Bethe Ansatz calculation 
and is 

e(p) = yAV- (5.44) 

We immediately notice that (15.441) differs from the free fermion case {epp{p) = 
^ p 2 ) just by a factor of A 2 and coincides with the latter (as expected) at A = 1. 

Therefore, we can repeat the exact same analysis as in the last section 
and derive the EFP (DFP) very easily. Introducing the complex Riemann 
invariants 

w = \ixp + iv (5.45) 
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and repeating the calculations of the previous section we obtain 

P DFP (R) ~ exp|-iA[7r(p -p) J R] 2 |, (5.46) 
P EFP (R) ~ exp|-^A(7rpoi?) 2 }. (5.47) 

Once again we can compare this result with what is already known ( TATBT) 
and confirming that, indeed, (15.47p gives the exact leading asymptotics of the 
EFP for the Calogero- Sutherland model. Subleading (in 1/R) corrections to 
(I5.47P are due to gradient corrections to the hydrodynamic action (I5.44p and 
to quantum fluctuations around the found instanton which go beyond our 
accuracy. 
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Chapter 6 

Hydrodynamics for a 
Spin-Charge System 

One of the most interesting predictions of the Luttinger Liquid (LL) model is 
the effect known as the spin charge separation, according to which, in a one- 
dimensional system, the fundamental excitations are "holons" and "spinous" 
and they carry independently the charge and spin degrees of freedom of the 
electrons respectively. 

These results of the LL model have been known for many years, but only 
now are we reaching the technological advancements necessary to test it exper- 
imentally. In recent years many laboratories have attempted to confirm this 
prediction (see, for instance [3]-[l3]) and the data are definitely in accordance 
with the theory. What is still missing is a conclusive test that would firmly 
link the experimental results to the theoretical prediction, in that the methods 
used so far are fairly indirect and their interpretation is not unique. 

While devices are developed right now to overcome this problem and confirm 
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with certainty the occurrence of spin-charge separation in one dimension, we 
are interested in looking ahead and investigating the correction to this predic- 
tion. 

In fact, it is well known that this effect is valid only for low energy excita- 
tions, in that its derivation assumes the spectrum to be linear. It is also known 
that the curvature of the spectrum will destroy perfect spin-charge separation 
for higher energy excitations and introduce a coupling between the spin and 
charge degrees of freedom. 

Unfortunately, the field theory description underlying the LL model (bosoniza- 
tion), becomes inconsistent if one were to consider a non-linear spectrum: even 
a reasonable quadratic spectrum generates a theory with no stable vacuum. 
Even if one tries to perturbatively include corrections to the linear spectrum 
approximation, the expectation values calculated with this theory are diver- 
gent and to complete a reasonable bosonization calculation beyond the linear 
spectrum approximation is still an open challenge in many contexts. 

To address the problem of determining the corrections to exact spin-charge 
separation we propose the hydrodynamic approach. We already discussed 
with the example of the EFP how the hydrodynamic approach is a natural 
generalization of a bosonization description to include the full spectrum of 
the microscopic theory. Therefore, we expect to be able to generalize the 
bosonization description of electrons in terms of non-interacting holons and 
spinons to include curvature corrections that will couple the two degrees of 
freedom. 

This hydrodynamic description will involve two interacting fluids, one de- 
scribing the bulk motion of the system, i.e. the charge, and one the internal 
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degrees of freedom, i.e. the spin. In the limit of linear spectrum the interac- 
tions between these two fluids will vanish and one would recover traditional 
LL results, while normally the interactions will mix spin and charge in the 
fluids. Ideally, we might find that the gradient-less hydrodynamic description 
is integrable in the sense of Appendix [Dl Then, one could be able to derive the 
stable excitations of the system, "solitons", i.e. one could identify two stable 
quasi-particles that would carry a fraction of the total charge and of the total 
spin of an electron each. 

To derive the hydrodynamic description, we use the exactly solvable model 
of electrons with contact interactions, which was first solved using the Bethe 
Ansatz in [50] • Another interesting integrable model to consider would be the 
spin Calogero-Sutherland model [51] , which can be solved using the asymptotic 
Bethe Ansatz. These two models are at opposite limits, in that one assumes 
the shortest interaction possible, i.e. contact interaction, while the other de- 
scribes an inverse quadratic long range potential. While these can be good 
approximations to some physical systems, most systems lie in between these 
two extremes. Nonetheless, studying the limiting cases, where the theory of 
integrable models can help us, will allow us to us gain some insights on the 
general structure of a hydrodynamic description of spinful electrons. 

In section 16.11 we sketch how one would bosonize a theory with quadratic 
spectrum and show that the resulting model has an unstable vacuum and that a 
perturbative treatment of the curvature of the spectrum generates divergences 
in the calculation of physical observables. In section l6~2| we describe the Bethe 
solution for a gas of electrons with contact interaction and we show how to 
implicitly derive its hydrodynamic description. 
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6.1 Bosonization for a quadratic spectrum 

When one wants to construct the bosonization description of a system, the 
first step is to linearize the spectrum. As we are going to show, it is hard to 
make sense of the bosonization of a non-linear spectrum. Therefore, often one 
refers to "bosonization" as the whole procedure of linearizing the spectrum and 
then expressing the system in terms of its density excitations. However, it is 
important to keep in mind that in this section we are going to distinguish these 
two steps and call "bosonization" just the transformation from the microscopic 
degrees of freedom to the collective ones. In doing so we are aware that we are 
losing rigorousness and that the resulting model has to be interpreted carefully. 

With these remarks in mind, let us now sketch how one would bosonize a 
quadratic theory. For simplicity, let us consider free fermions: 

U = -¥(x)d 2 x ^(x) 

= k 2 m\k)m{k) (6.1) 

where d x = d/dx and the second line shows the spectrum of the Hamiltonian 
in the Fourier space representation. 

The traditional bosonization approach linearizes the spectrum as 

L,R 

~ - k FY^L,R^ L ' RTikF J2^L,R d ^L,R + ■■■ (6.2) 
L,R L,R 

where the first term is interpreted as a chemical potential and can be absorbed 
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in a redefinition of the ground energy, while the second term expresses the 
linear spectrum of the excitations around the Fermi points ±kp. Moreover, 
we defined left- and right-moving fields ipL,R as the fields obtained expanding 
around the left /right Fermi Point: 

V(x) = : e ikpX i; R (x) + e-' lkFX ^ L {x) : . (6.3) 



In bosonization one effectively describes the system with a single collective 
field that captures the density fluctuation. In practice, one makes the following 
transformation: 



1 



27T 



(6.4) 



where 4>l,r are the left and right moving bosonic fields. 

We use the transformation ( 16.41) to calculate various bilinears in the spin 
fields. For instance, one can consider a quantity like 



_ . e ±iV^{<t>L,R{x+t)-4>L,R,{ x )) ■ e ±4-K{<f>L,R{ x )<i>L,R{ x +£)} 

1 



2i7re 



^±iV4n(<f>L,R{ x +£)-<l>L,R{ x )) _ J 



(6.5) 



where we used the identity 



e A e B = : e A+B : e ( AB + A ^) 



(6.6) 



and the fact that 

(^ji(O)^(x) - ^(0)> = lim J- In (6.7) 
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The colons denote normal ordering^ and in the second line of (16. 5p we used the 
fact the normal ordering simply amounts to subtract 1 from the exponential. 
We can expand (16.51) in powers of e to find 



n=Q U - 

= JL r ±iV5F££=i _ x ] (6S) 

2ivre L J v ' ; 

which give the generating function of the fermionic currents 

jt R (x) ee ^ LjR (x)d2^(x) (6.9) 

in terms of the bosonic field 4>l,r- 

By matching powers of e in H 6 . 8 f) we can write down these expressions. The 
density of fermion is 



PL,R = jo ,R = i'[ R (x)lp L ,Ii( X ) = ±-^d x <f) L>R {x), (6.10) 



the current density is 



h' R = ^l, R (x)d x ^ L A x ) = i (d x (p LtR (x)) 2 ± —==^> LR . (6.11) 

V47T 



The third term in the expansion can be identify with the original quadratic 



1 Normal ordering means that in the evaluation we should put the creation operators 
to the right of all the annihilation operators. This convention is equivalent to subtracting 
un-physical zero energy contributions. 
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Hamiltonian for the left /right movers 



(6.12) 



but is important to notice that while the first line is a well defined Hamiltonian 
operator, the second line is not and can be understood only as a perturbative 
interaction term. 

The last terms in (16.121) are total derivatives and can therefore be neglected 
as boundary terms. Therefore a model of free fermions with quadratic spec- 
trum would be transformed into a theory of bosonic fields with Hamiltonian 



This describes a cubic theory and therefore it cannot be quantized, since the 
spectrum for the bosonic field has no lower bound and the ground state of the 
theory is unstable and has an infinite energy. This is the reason for which it 
does not make sense to directly bosonize a non-linear theory. 

If one were, instead, to consider the linearized version of the fermion theory 
( 16. 21) . using the expressions found above, the bosonized Hamiltonian would be 



One can then include the additional terms like (16. 12ft neglected in (16. 2h in 



n = (d x <f )R f + (d x <p L f. 



(6.13) 



H ~k F (d x (j> R f + k F (d x (f) L ) 2 + ... 



(6.14) 



122 



a perturbative way and treat them as small correction. Unfortunately, even 
this attempt is ill-fated, since the calculations of observable quantities diverge 
and nobody has found the correct way to resum the diagrams and cure these 
infinities. 

To add spin degrees of freedom in the bosonization is very easy. Essentially, 
the machinery we just outlined can be repeated for two fields, corresponding to 
spin-up and spin-down fermions/bosons. The fermionic version of the Hamil- 
tonian will therefore be 

H = M' (.r)r7;M' (./•) - *\(x)%^(x) (6.15) 

while the linearized bosonic one will look like 

H^k F j2(9A, R y+kFj2{ d ^U) 2 + --- ■ ( 6 - 16 ) 

L,R L,R 

One of the most interesting properties of the bosonization procedure is 
that all the interactions between fermions that are relevant, i.e. that cannot 
renormalized to zero, can be expressed in terms of quadratic interactions for 
the bosonic field. One can show that the effect of the interactions can be 
absorbed in a different value for the coefficients in front of the left and right 
moving fields in (16.141) . 

One can define spin and charge bosonic fields as 

C = T + 0^ 

S = T -0^ (6.17) 
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which are often referred to as "holons" and "spinous" respectively. In terms 
of these fields, the Hamiltonian is 

H ~ v c (d x <j) c ) 2 + v s (d x <j) s ) 2 + ..., (6.18) 

where we took into account that interaction can renormalize the coefficients 
in front of the spin and charge degrees of freedom differently, redefining the 
spin and charge (Fermi) velocity. 

The Hamiltonian being quadratic, transformation (16.171) leaves it substan- 
tially unchanged, i.e. (I6.18P remains quadratic. But if we were to include the 
additional cubic terms coming from the spectrum curvature, transformation 
(I6.17P would mix the spin and charge degrees of freedom, effectively introduc- 
ing interactions between holons and spinons like 

AH = ^- (d x <J) c ) {d x <t> s ? + j- {d x <t> c f (d x cf> s ) . (6.19) 

Kp Kp 

These terms are suppressed compared to the main ones in (16 . 1 8[) by a factor 
of 1/kp (throughout the calculation, kp has been the "big" parameter in the 
theory). One could therefore introduce them in the bosonization description 
and use them to calculate the corrections to exact spin-charge separation. 
Unfortunately, such calculations are again ill-fated in that their results are 
divergent. This is quite normal in bosonization calculations, but, as far as we 
know, nobody has been able to devise a scheme to cure this divergences for 
this problemc- 

2 Let us remark that recently Pustilnik and coauthors have successfully solved a similar 
problem in [52] . 
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That is the reason for which we think that a hydrodynamic approach to 
the problem can be quite effective in finding the corrections to spin-charge 
separation at low energies. 



6.2 Hydrodynamics for fermions with contact 
interaction 

We now consider a system of N spin-1/2 fermions interacting through the 
Hamiltonian 



where c is the coupling determining the strength of the contact interaction. 

We want to consider the Bethe ansatz ground wavefunction for the state 
with a fixed number of spin down particles M, with a fixed total momentum P 
and a fixed total spin momentum P s (defined, in this the momentum 

of the spin-down degree of freedom). In a Grand Canonical approach we would 
write: 



H = H + v h {M) + A t 1 / il (2M -N) + fi 2 h 2 (P) + V 3 h(P s ). (6.21) 



These definitions are chosen so that N and P, and M and P s are canonically 
conjugate hydrodynamics variables satisfying a continuity equation for each 



The Bethe Ansatz construction is more complicated than the one outlined 
in Appendix O because in this case we have two species of particles coexisting, 




(6.20) 



pair. 
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fermions with spin up and fermions with spin down. Therefore, one needs 
to include an additional degree of freedom in the Bethe Ansatz approach, a 
degree of freedom corresponding to the spin of the particle that is going to be 
parameterized by a spin quasi-momentum, or, to be more precise, by a spin 
rapidity^. 

The existence of this second sets of parameters for the spins brings a second 
set of Bethe equations ( 1C.14|) that have to be satisfied. In the thermodynamic 
limit this will produce two coupled integral equations in the density of the 
quasi-momenta and in the density of the spin rapidities. The Bethe Ansatz 
construction for this model was firstly reported in [50], while its thermody- 
namic study was undertaken in [53] . 

We assume a wavefunction of the form 

*(xisi, x 2 s 2 , . . . , x N s N ) = $f Oi, x 2 , ■ ■ ■ , x N )Gf (6.22) 

j 

where = ±1/2 are the spin quantum numbers and Gj 1 is the spin part of 
the wavefunction. The spatial part of the wavefunction is constructed as 

§f = [Q, P]e l E ^ k ™ XQ > (6.23) 
p 

with the coefficients given by 

M 

[Q, p] = ± Yl A * II y*) ( 6 - 24 ) 

R 3=1 



3 Rapidities are variables related to the quasi-momentum that allow to express the scat- 
tering matrix in a simpler form. 
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where 



kpn — A + ic 



fJl k PU+i) - A - ic 
An = TT ! flJ "! K + 2iC - (6.26) 



The ?/j's are the coordinates of the spin down particles (i.e. a subsets of the 
Xj's); Q, -P and i? are respectively permutations of the Xj's, fc/s and t/j's. 

By imposing periodic boundary conditions we find the following Bethe 
equations 

Ka — A a + ic 



e i kjL = tt y ^± (6 .27) 



0=1 J 

jV 



n A a - fcj + ic _ rr A« - A/3 + 2ic 
^-^•-ic " 11 A. - A, - 2ic (6 ' 28) 

constituting a nested Bethe ansatz. Taking the logarithm of these equations 
we obtain 

Lk = 2nl k - 2 arctan — ^ (6.29) 

= 2ttJ a - 2 ^ arctan ( ^— - j + 2 ^ arctan ( A - A J . (6.30) 
A V c / A , V 2c J 

Due to the additional degrees of freedom, we now have two sets of integers Ik 
and J\ to define the state of the system. 

In the thermodynamic limit (L — > 00, N ^ 00, M — > 00 with the condition 
that p = N/L and p s = M/L are finite) the sums can be converted into 
integrals and the densities of momenta fc's and rapidities A's are determined 
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through integral equations: 

[ Br 4ccr(A')dA' [ Qr 2cr(k)dk 
f BR 4ca(A)dA 

2 ^> = 1+ J BL c- + V-Ar (6 - 32) 

These two integral equations define the densities r(k) and cr(A) as a func- 
tion of the parameters Ql,Qr, Bl, Br. To determine these parameters in 
terms of physical observables we need to satisfy the following consistency con- 
ditions: 



p= r(k)dk, P= kr(k)dk, 

JQl JQl 

Ps = / cr(A)dA, P s = / p(A)<r(A)dA (6.33) 
Jb l Jb l 



where p 
pidity A|f 



A) is the pseudo-momentum of the spin degrees of freedom with ra- 



<-Qr f\ r . _i_ a _ u s 
p(A) = iM / In --— ) r{k)dk , (6.34) 



L 



ic - A + k 



which comes from the definition 



N . , . 

7=1 J 



Finally, the energy (Hamiltonian) of the system is given by 

H(Q L ,Q R ,B L ,B R )= I k 2 r{k)dk. (6.36) 

Ql 



4 These definitions give P = ^ J2k + TT and P s = 2£ ^ A J A in terms of the 

Bethe integers. 
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We can now invert the system of equations (16.331) to express the parameters 
Qli Qr, Bl, Br in terms of the hydro dynamic variables p, J, p s , J s and to con- 
struct the hydrodynamic Hamiltonian 

H(p,P,p s ,P s ) = H(Q L ,Q R ,B L ,B R ). (6.37) 

The equations of motion can be easily derived using the fundamental commu- 
tation relations 

\pi(x),vi(v)] = IPlWMv)] = -&{x-y) (6.38) 

where v = j/p is the velocity and the spin/charge degrees degrees of freedom 
are connected to the spin-up/down one by 

p = P] + Pl J=J^ + J l (6.39) 

Ps = Pi Js = Ji- (6.40) 
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Chapter 7 

Aharonov-Bohm effect with 
many vortices 

We now turn to a problem that is very different from what we considered 
so far and we study a two-dimensional configuration. The Aharonov-Bohm 
effect is the prime example of a zero field situation where a non-trivial vector 
potential acquires physical significance, a typical quantum mechanical effect. 
We consider an extension of the traditional A-B problem, by studying a two- 
dimensional medium filled with many point-like vortices. Systems like this 
might be present within a Type II superconducting layer in the presence of a 
strong magnetic field perpendicular to the layer. We are going to construct an 
explicit solution for the wavefunction of a scalar particle moving within one 
such layer when the vortices occupy the sites of a square lattice. From this 
construction we infer some general characteristics of the spectrum and imply 
that such a flux array produces a repulsive barrier to an incident low-energy 
charged particle, so that the penetration probability decays exponentially. 
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7.1 Introduction 



In classical mechanics it is said that the vector potential has no physical mean- 
ing, and only the electromagnetic field has physical (measurable) effects, since 
it is the only gauge invariant. In quantum mechanics, however, the vector 
potential appears in gauge invariant quantities that describe a new class of ef- 
fects. In these cases, corresponding to topologically non-trivial configurations, 
we recognize the importance of the vector potential, even when the electromag- 
netic field vanishes everywhere in the regions accessible to a charged particle. 

The standard example of this class of effects is recognized in the Aharonov- 
Bohm effect [51], in which a magnetic field is confined to a region of space, 
and electrically charged particles are only free to move outside this region. 
Although a particle cannot experience the field strength directly, the covariant 
momentum, i.e. the momentum derived through the minimal coupling to the 
electromagnetic field, is 



and is affected by this overall configuration, because the vector potential 
carries a 'memory ' of the presence of the magnetic field even outside the region 
where the field is localized. 

In this way, the particle is influenced by the field, through a shift in the 
phase of the wavefunction 



where $ is the total magnetic flux inside the circuit (i.e. a closed path of the 





(7.2) 
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particle J*. This explains why the effect is called 'topological': the behavior of 
the particle is sensitive to the overall configuration of the system, even though 
there is no classical magnetic force at any point. 

The extension of the A-B problem in the presence of many localized fluxes 
cannot in general be tackled exactly. There exists a simple argument [56] 
due to Aharonov which shows, using the Bloch theorem, that an infinite line 
of equispaced point-like fluxes would constitute an impenetrable barrier to a 
particle of sufficiently low energy. The particle would not be able to pass 
through such an array because it could not satisfy simultaneously on both 
sides of the barrier the Bloch periodicity conditions on its phase, in the light 
of the A-B effect. 

We are interested in exploring a possibly more realistic set-up by studying 
the propagation of a charged particle through a medium filled with point-like 
fluxes. 

Experimentally, one might find a situation similar to this inside a Type II 
superconducting layer in the presence of a large magnetic field perpendicular 
to the layer. Quasiparticles in the layer would encounter numerous vortices, 
each containing a superconductor flux quantum, and under some conditions 
might not penetrate the vortices (see, for instance, [57]). 

Situations similar to this have been addressed by several authors in recent 
years [58], especially in connection with transport studies and with a special 
focus on the Hall conductivity of two-dimensional electron gases on top of 
superconducting material. These works have shown a depletion of the density 



^or a phase shift to be measurable, we need to create an interference effect and therefore 
the particle has to come back to the same point. It does not make sense to speak of a phase 
shift for open paths. 
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of states at the bottom of the spectrum. 

Our aim is to consider such a 2-dimensional layer, punctured by magnetic 
fluxes, and to study the wavefunction of a single scalar particle entering this 
medium. For simplicity, we take the vortices as point-like, so that the space 
available for the particles is a punctured plane. A similar attempt was done 
by Nambu in [55], and our aim is to extend his results. 

We are going to show that a lattice of impenetrable magnetic fluxes (vor- 
tices), such as the one described above, constitutes a barrier to a low-energy 
charged particle trying to pass through the medium. That is, the distribution 
of the vortices creates a configuration whose topological constraints on the 
wavefunction are comparable to an effective repulsive potential. Qualitatively, 
there number of ways to see this: 

• The presence of the fluxes generates a non-zero vector potential inside the 
medium, raising the minimum energy (that is the square of the covariant 
momentum, Eq. 17. 2p required for an electrically charged particle to exist 
in the medium, 

• Particles are repelled by the vortices, as their wavefunctions must vanish 
on the vortex sites. Therefore, the bigger the typical amplitude of the 
wavefunction in the flux-containing region, the bigger the energy due to 
the sharp spatial variation. This means that for low-energy states the 
wavefunction will not be able to reach a value appreciably different from 
zero in the presence of fluxes, 

• The analysis of Nambu indicates that the medium constitutes a barrier 
even from the point of view of angular momentum. In [55] , he argues that 
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the angular momentum of a particle should be greater than the magnetic 
flux present in the medium if the particle wavefunction is to satisfy the 
boundary conditions. In other words, the lower angular momentum levels 
are missing and are not part of the spectrum. 

This result is even more clear in light of the aforementioned works [58] 
showing a Lifshits tail in the density of state at low energies for a random 
distribution of vortices. From a physical point of view, it seems quite clear 
that a charged particle approaching the medium with sufficiently low energy 
will be repelled, that is, its penetration will be exponentially damped. In 
the same way, if we localize a particle in its ground state in a region without 
vortices, the particle will not be able to escape outside that region through one 
containing vortices except by tunneling, and we should be able to construct 
a bound state of topological character (actually a very long-lived resonance), 
even though there is no classical force. The fact that a bound state can be 
topological in nature is new and was already suggested by the work of Nambu 



7.2 Mathematical preliminaries 

We concentrate on the case in which all the N fluxes have equal strength 
$ = $o/2, where $0 = 27r| is the quantum unit of flux. In this case it can 
be shown (see, for instance, [59]) that the problem is invariant under time- 
reversal, and we can therefore choose the wavefunctions to be real. 

Indicating with (xi,yi), i = 1..N, the coordinates of the vortices, we can 
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write the vector potential in the standard circular gauge as 



(N N 
v ^ y — Hi \^ x — Xi 

( X _ x .)2 + {y _ y .)2 ' - (x _ x . )2 + (y _ yi) 2^ 

N 



i=i 

N 



3=1 

N 



vxA = 2tt$ 62 ( x ~ x i^y ~Vi)- ( 7 - 4 ) 



i=l 



The equation of motion for a particle in this medium is given by the 
Schrodinger equation (in units h = e = 1) 



— (V-iA) 2 * + £* = 0, (7.5) 
2m 



and, in these units, integer $'s correspond to a quantized flux (in our case 
$ = 1/2). 

Following the idea of Nambu [55] , we implement a singular gauge transfor- 
mation G to remove the vector potential: 



c -n( i::g;MS )-- 



In this way, we reduce our problem to a free-field case 



- -^V 2 ^ = Erf,, (7.7) 
2m 



with non-trivial (topological) boundary conditions on the wavefunctions in the 
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region surrounding each vortex. 

In constructing our solutions, we must require that the wavefunctions van- 
ish on the vortex sites 

^(x = Xi , y = Vi) = i = l..N , (7.8) 

and that they acquire the A-B phase e 2l7r * = — 1 each time a particle completes 
a turn around a vortex. More precisely stated, in this singular gauge the effect 
of the vector potential is represented by a phase-matching condition on the 
wavef unction 

^(0) = -i/j(9 + 2tt) (7.9) 

where 9 is the azimuthal angle about the vortex. 

We know from standard complex analysis that this condition implies the 
existence in the 2-dimensional plane of a cut connecting two distinguished 
Riemann sheets. For a real wavefunction this last condition implies that there 
exists at least one line exiting each vortex site on which the function has to 
vanish in order to change its sign. 

7.3 Construction of the solutions 

The construction of the solution on a general distribution of fluxes is not easily 
attainable (as argued in [55]). We do not need to confront these complications 
in order to show our point and so we shall simplify the problem by taking 
the vortices as located on the vertices of a square lattice of lattice spacing L 
(Fig. 17.11) . a case for which we shall be able to give an explicit solution to the 
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Figure 7.1: The vortices are located 
on the sites of a square lattice. 



Figure 7.2: The vortices are paired 
and connected by segments on which 
the wavefunction has to vanish in or- 
der to satisfy the topological condi- 
tions. Grey lines indicate the real pe- 
riodicity of the lattice and identify the 
fundamental region over which we will 
work. 
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problem. 

Inspired by a recent construction [60], we try to give an estimate of the 
minimal energy required for a charged particle to exist in the medium, and 
also to calculate the decay factor of particles with zero energy in the lattice. 

Before we construct the solution, it may be helpful to spend a few more 
words on our boundary conditions. Since we can take the wavefunction to be 
real, we translated its phase shift around each vortex with the condition that 
the solution has to vanish along one line, but we have not specified this line. 
This line is not the familiar cut in a complex plane (which is, of course, a 
gauge choice). In fact, we have some freedom in the choice of the line along 
which the wavefunction vanishes, but this is not a gauge freedom in that it 
has a measurable effect. It would be better to say that the position of this 
line is a freedom of choice for the wavefunction. Therefore, in order to impose 
it as a boundary condition, we have to make this choice appropriately for the 
problem we want to study (this consideration will be important when we will 
consider the penetration of a zero-energy solution inside the medium). 

Let us consider for a moment just a pair of vortices. If we choose the 
line on which the wavefunction has to vanish as the ray exiting one vortex 
and pointing in the direction of the other one, we can see that the boundary 
conditions become that the function has to vanish only along the segment 
connecting the two fluxes; this is certainly a very convenient choice, compared 
to other solutions which would require the wavefunction to vanish on two 
semi-infinite lines and therefore to develop higher gradients. 

To construct the lowest energy solutions let us consider the vortices in pairs, 
connecting nearest neighbors with line segments along which the solution has 
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Figure 7.3: Boundary conditions and 
parity of the wavefunction: the black 
continuous lines represent Dirichlet 
boundary conditions, while the grey 
dashed lines indicate Neumann condi- 
tions. 
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Figure 7.4: The region over which we 
construct the fundamental solution. 
The rest of the lattice can be covered 
starting from this basic tile. The con- 
tinuous black line indicates where the 
wavefunction must vanish (Dirichlet 
condition) and the dashed ones where 
its derivative is zero (Neumann con- 
dition). We expand the solution on a 
basis in the region I and on another 
basis in the region II and we impose 
continuity of the function and deriva- 
tive across the grey line. 
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to vanish. For definiteness, we connect fluxes on the horizontal direction, 
requiring the wavefunction to change sign when it crosses these segments (see 

Fig. pp. 

Along these segments the wavefunction possesses odd parity. If we are 
interested in the low energy modes, this means that along the continuation 
of these segments, the function will be even and so its derivative must vanish 
there. To conclude our analysis on the boundary conditions, we notice that 
our system is clearly periodic. To ensure periodicity of the wavefunction, 
we require its derivative to vanish identically along the sides of each square 
centered on a flux (see Fig. 17. 31) . 

Bearing these considerations in mind, we now have to solve a problem with 
mixed Dirichlet and Neumann boundary conditions. We can further reduce 
the system under study and concentrate on two of the quadrants around a flux 
site, because the rest of the lattice can be covered by mirroring and flipping 
this unit (Fig. 17.41 in reference to Fig. 17.31) . 

In summary, we now have to solve the problem of a free particle in a 
rectangular box with sides of length 2 and 1 (in units of half of a lattice 
spacing). We impose Neumann boundary conditions everywhere, except on 
half of one of the long sides, where we require the Dirichlet boundary condition. 

This is a non-standard problem; as we are not aware of any previous study 
on a system with these boundary conditions, we shall proceed in constructing 
the solution starting from a basis compatible with the conditions. In region I 
of Fig. 17.41 we identify a convenient basis in the set 

{cosh [k n (l + x)] cos(niry)}™ =0 , (7.10) 
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Figure 7.5: We truncate the infinite-dimensional matrix to a size N and we find 
the first energy eigenvalue Eq = 2mEo corresponding to this shorten system. 
This is the plot of N versus Eq and its fit with a polynomial in inverse powers 
of iV up to the third order (higher orders do not contribute appreciably). 
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while in region II we expand the solution on 



oo 



cosh [i£n(l — x)] sin (n H — )iry 



(7.11) 



n=0 



with the condition n 2 ir 2 — k\ = [n + |) 2 vr 2 — K 2 = 2mE. 

By matching the wavefunction and its derivative across the line x = 0, we 
may seek the values of e — 2mE for which the system admits a solution. In 
principle, this would involve the calculation of the determinant of an infinite 
matrix. To obtain an approximate solution, we truncated the system to a 
finite size, and found the first energy eigenvalue Eq = 2mE as a function of 
the size of the matrix (see Fig. 17.51) . Then, we plotted Eq versus the order N 
of the matrix and performed a fit with a polynomial in inverse powers of N, 
taking the zeroth-order coefficient as the solution we would have obtained by 
considering the whole infinite system. 

In this way, we found the first energy eigenvalue for our solution to be: 



Next, we are interested in estimating the decay factor of a particle enter- 
ing the medium with zero energy. This problem depends on the direction in 
which the particle is traveling, in that it is connected with the choice of the 
ray/segment over which the solution has to vanish. Depending on the direc- 
tion of motion, the wavefunction may 'choose' different configurations for these 




(7.12) 



that is 



E = (2.0682 ± 0.0002)m" 1 L- 2 



(7.13) 
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Figure 7.6: Decay of the zero-energy solution moving horizontally. We re- 
quire periodicity on the vertical axis and exponential decay in the horizontal 
direction. The continuous black line indicates where the wavefunction must 
vanish (Dirichlet condition) and the dashed ones where its derivative is zero 
(Neumann condition). 



We solve the problem for a particle moving along the x direction. That 
is, we construct a solution which exhibits periodic behavior in the y direction 
and real decay in the x (Fig. \7.6h . Again, we expand the wavefunction on 
appropriate bases: in region I and III of Figure 17.61 we use 



segments. 



{e 



nnx 



cos(ran/)} 



(7.14) 



n=0 



for right-moving and 



{e 



nnx 




oo 



(7.15) 



n=0 



for left-moving modes. In region II we expand on 




(7.16) 
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Figure 7.7: We truncate the infinite-dimensional matrix to a size N and we 
find the lowest value for the decay factor K corresponding to this finite system. 
This is the plot of N versus K and its fit with a polynomial in inverse powers 
of iV up to the third order (higher orders do not contribute appreciably). 



for right-moving and 



sin 



in 



(7.17) 



n=0 



for left-moving modes. 

We impose matching of the wavefunction and its derivative across the line 
x — — 1 and x = 1 and we write the damping of the solution by requiring an 
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exponential suppression: 

^(x = -2,y) = e 4 ^(x = 2,y) ^-(x = -2,y) = e 4K ^(x = 2,y). 

(7.18) 

We look for the values of K for which the system admits solution. 

As before, the system of equations is infinite-dimensional, so we found the 
lowest value for K as a function of the order N of the matrix and performed 
a fit with inverse powers of N to retain the zeroth order of the polynomial as 
the solution (see Fig. 17.7ft . 

In this way, we find a decay factor for a particle moving along the horizontal 
direction: 

K = (0.88 ± 0.01) x - = (1.76 ± 0.02)L" 1 , (7.19) 

and the same K holds for a particle moving in the vertical direction because 
we have the freedom to rotate the system by 90 degrees and rearrange the 
segments connecting the vortices in the new direction. 



7.4 Conclusions 

Considering a lattice of point-like magnetic vortices, we showed that the spec- 
trum for a particle in such a medium is discrete and that the lowest energy 
eigenvalue is greater than zero, by explicitly constructing the wavefunction 
(E = (2.0682 ± 0.0002)m~ 1 L- 2 ). 

This is quite in contrast with what was predicted by Y. Nambu in [55]. In 
this paper, the author argues that a solution of the Schrodinger equation in 
our gauge would have to be either holomorphic, or anti-holomorphic. 
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To see this, let us switch to complex coordinates to describe the plane. The 
free particle equation now reads: 

d t d 3 il> = EiP (7.20) 

and therefore the solution for zero energy is either analytical or anti-analytical. 
Nambu argues that, by continuity, this property should persist at higher ener- 
gies as well. In the preceding section we constructed a non-zero energy solution 
which clearly is neither holomorphic nor anti-holomorphic. 

The analyticity or anti-analyticity of the solutions is an important point 
of Nambu's construction that leads him eventually to conclude that the states 
with lower angular momentum are not admissible in the spectrum. This would 
imply that a particle entering the medium with zero energy would undergo a 
suppression, which is not merely exponential, but at least Gaussian. For that 
reason, we argue that our approximation comes closer to the true behavior, 
because by allowing more penetration it reduces uncertainty-principle energy. 
This statement applies even for zero energy, where Nambu's argument appears 
rigorous at first sight from (I7.20p . The loophole, we believe, is that in any case 
the wavefunction is not analytic at the location of a vortex, and therefore the 
factorization of the Laplace operator fails at that set of points. 

We computed the decay factor for a particle moving along one of the lat- 
tice directions to be K = (1.76 ± 0.02)L~ 1 , and showed that this decay is 
purely exponential. The magnitude of this suppression depends on the direc- 
tion of travel. To compute the decay factor in other directions it would be 
necessary to modify ad hoc the boundary condition (the positioning of the 
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ray where the wavefunction vanishes). The condition we worked with is the 
one that minimizes the extension of such rays and therefore poses the minimal 
constraint on the solution. Any other choice would have a greater impact on 
the shape of the wavefunction, and this suggests it would shorten the decay 
length. The directional dependence is easy to understand, because the cou- 
pling between charge and vortex is strong, so that the lattice length scale and 
the decay length are comparable: in the limit of vanishing lattice constant 
the decay length also vanishes. A quantitative analysis for generic directions 
would require a different formalism from the one implemented here. 

One can of course estimate the lowest allowed energy level with a mean 
field approximation. A long wavelength particle (k ~ 0) would see a uniform 
magnetic field 

and it would occupy the lowest Landau level 



E— = hiu c - = — (7.22) 
2 mc 2L 2 v ; 



This naive calculation disagrees somewhat with the value we found for the 
lowest energy level and indicates that is probably possible to construct a better 
wavefunction: 

Tp rp naive 

° * = 24% (7.23) 

Nonetheless, we believe that the order of magnitude of the effect has been 
established, in that the energy eigenvalue and the decay rate K agree in this 
respect. This agreement implies that the topological constraints imposed by 
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the configuration of vortices act as an effective repulsive potential of order 
unity (taking the mass m — 1), and that this potential is direction-dependent. 
This effective potential could be used to trap a particle in a region, just by 
surrounding that region with a medium of localized fluxes. This may be a new 
form of trapping. 

It would be interesting to investigate the more physical situation in which 
the vortices are not forced to occupy the sites of a square lattice, instead being 
randomly distributed, as in [58] and maybe dynamical objects themselves, 
but this would require more powerful mathematical tools. It seems plausible 
to us that the order of magnitude of the decay length and the qualitative 
characteristics of the problem would not be very different from the ones found 
with our model. Our reason for saying this is that one could replace the 
random vortex distribution with a random distribution of short line segments 
on which the wavefunction vanishes, and this array surely would be equivalent 
to a repulsive potential of characteristic magnitude, leading to exponential, 
not Gaussian decay. 
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Chapter 8 
Conclusions 



During my Ph.D. studies, we have devoted considerable efforts to the study of 
the correlator known as "Emptiness Formation Probability" (EFP). This cor- 
relator measures the probability that a one-dimensional system spontaneously 
develops a region of space depleted of particles. 

This correlator is a n-point correlator, where n is the number of particles 
that one has to move to empty the region of space, so its complexity grows with 
n. Despite this apparent complexity, the EFP was introduced as a very natural 
correlation function in the development of a determinant representation for 
correlators of integrable models using the Bethe Ansatz technique. 

The calculation of correlation functions for exactly solvable models is an 
arduous task and while some exact expressions have been derived in terms 
of operator-valued determinants, they are often of very limited practical use. 
The EFP has the simplest expression in this representation and is therefore 
believed to be the 'fundamental' correlator for integrable theories. 

In Chapter [3] we calculated the asymptotic behavior for n — > oo of the 
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EFP in the anisotropic spin 1/2 XY model. We showed that for this system 
the EFP has an exact expression in terms of the determinant of a complex 
matrix belonging to a family of matrices known as "Toeplitz Matrices". We 
used some results known in the theory of Toeplitz determinants to calculate 
the asymptotic behavior of the EFP throughout the rich phase diagram of the 
XY model. We found that the EFP decays exponentially almost everywhere, 
except on the isotropic line, where the decay is Gaussian. We also identified a 
power-law prefactor with universal exponents on the critical lines of the phase- 
diagram. Using a bosonization approach, We interpreted the crossover from 
Gaussian to exponential behavior as an intermediate asymptotics effect. 

Beside its importance in the theory of integrable models, the EFP is an 
interesting correlation function in that it measures a substantial deviation 
from the equilibrium state of the system that requires contributions even far 
from the Fermi points. For this reason, the EFP is a perfect test field for a 
hydrodynamic description of integrable theory. 

Very often, calculations of correlators in one- dimensional integrable sys- 
tems are performed using the bosonization technique. While being very pow- 
erful, the bosonization approach has a strong limitation in the assumption of a 
linear spectrum. Because of this, calculations done within this approximation 
are valid only as long as only excitations close to the Fermi points are involved. 
This assumption is clearly violated by the EFP. 

We introduced the hydrodynamic description of integrable systems and we 
showed its ability to correctly calculate the leading asymptotic behavior of the 
EFP for a couple of Galilean Invariant systems (Free Fermions and Calogero- 
Sutherland particles), as firstly derived in [19] . 
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In [19], the hydrodynamic description for free fermions on a lattice was also 
considered and its EFP was calculated to leading order. We have not included 
this project in this thesis, but the problem of applying the hydrodynamic 
approach to other lattice systems is definitely at hand. One should be able to 
derive the EFP of the XY model and of the XXZ model from a hydrodynamic 
point of view. 

In Chapter [6] we addressed the prediction of exact spin-charge separation 
coming from the Luttinger Liquid model. We argued that this result is a 
feature of the linear spectrum approximation and that any curvature of the 
spectrum would lead to a coupling between charge and spin degrees of freedom. 
We also explained some difficulties encountered by the bosonization method 
to account for these corrections. 

A hydrodynamic description naturally takes into account the interactions 
between spin and charge degrees of freedom. We sketched the derivation of the 
hydrodynamic description of fermions with contact repulsion. The hydrody- 
namic description is given in an implicit form, leaving the analysis of the cou- 
pling between the spin and charge degrees of freedom for further work. Another 
interesting model to study in this respect is the spin Calogero-Sutherland. 
Fermions with contact repulsion and this latter model represent the two inte- 
grable extremes of short and long range interaction, respectively, and therefore 
it is important to understand how spin and charge couple in both systems. 

Another open question is the origin of the double infinite series of con- 
served quantities for gradient-less, Galilean invariant, hydrodynamic systems, 
as discussed in Appendix [D] The existence of these conserved currents was 
discovered in the 1980's and we believe to be the first to actually construct 
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them explicitly. It is still unclear what underlying symmetry is responsible for 
these conservation laws. 

Finally, in Chapter [7] we looked at a two-dimensional system in which a 
scalar particle moves into a medium filled with point-like magnetic vortices 
pinned at the sites of a square lattice. We studied the Aharonov-Bohm prob- 
lem connected with such a configuration when the magnetic fluxes are all equal 
to a half of the quantum flux unit. We explicitly constructed a wavefunction 
for this system, showing that the spectrum is discrete, and considered the de- 
cay of a zero-energy particle within this medium and showed that this decay is 
exponential. We proposed that this medium acts effectively as a repulsive po- 
tential of topological nature (since no real field is present outside the vortices) 
and that one could use such a configuration to confine (trap) electrons. This 
is a novel idea and it would be interesting to extend this analysis to electrons 
with spin, since they would couple to the vortices differently depending on 
their spins. 



152 



Appendix A 

Exact results for EFP in some 
integrable models 

In this appendix we are going to recapitulate some known results obtained for 
the Emptiness Formation Probability P{R) in various integrable one-dimensional 
systems as they were reported in [19]. We present the results as 

S = -\nP{R). (A.l) 

for brevity, but also to facilitate the comparison with the instanton action S op t 
we are going to introduce in the main text. In this appendix we always set 
m — h—1. 
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A.l Free continuous fermions 



For Free Fermions we can use the fact that the ground state wavefunction 
(more precisely |^/| 2 ) coincides with the joint eigenvalue distribution of the 
unitary random ensemble. For the latter the probability of having no eigen- 
values in the range 2R of the spectrum was obtained in [36] (see also |61j). 
The first few terms of the expansion in 1/R are 

S = \s 2 + \ In a - (± In 2 + 3C'(-1)) + 0(^ 2 ), (A.2) 

where we introduced the parameter 

s = irpoR. (A.3) 

A.2 Calogero-Sutherland model 

The Calogero-Sutherland model [HI 09] with A^-particles is defined as 



where pj = —id/dxj is the momentum operator of j-th particle and A is a 
dimensionless coupling constant (we will come back to the Calogero-Sutherland 
model in section [531) . For A = 1 we have free fermions, at A = free bosons, 
while in the case of general A the model (1A.4|) describes interacting particles 
with fractional statistics. 
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The ground state wavefunction is 



*G = l[(x j -x k )\ (A.5) 

j<k 

Using the form of the ground state wavefunction and thermodynamic argu- 
ments EH one obtains 



S = £ (np R) 2 + (1 - X)np R + 0(lni2). (A.6) 



or, by defining 



and 



we have 



s = VXwpoR (A. 7) 



A 1 / 2 -A- 1 / 2 
a = , (A.J 



S = \s 2 - 2a s + O(lns). (A.9) 



The notation ao originates from conformal field theory, since the theory 
with central charge c = 1 — 24a 2 , is known to be related to the Calogero- 
Sutherland model [62] • We are not aware of a determination of the coefficient 



in front of the In s term in the expansion for general A. However, A = 1/2, 1, 2 
correspond to random matrix ensembles where the full asymptotic expansion 
is known (see below). In those cases the coefficient of Ins is 1/8,1/4,1/8 
respectively. The natural guess for an interpolation is [19] 



S = -s 2 -2a s+ \ --al ) Ins + 0(1). (A.IO) 



2 V4 
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A. 3 Random matrices 



For Random Matrix ensembles with (3 = 2A = 1,2,4 the joint eigenvalue 
distribution is proportional to Yli<j \ z i ~ z jf ■ The ^ un asymptotic expansion 
of the quantity -E^O, 2R) corresponding to the EFP P(R) was obtained using 
properties of Toeplitz determinants (36j EI]. The first few terms of these 
expansions are given by 



Sx=i/2 = ^ 2 + ^ + lln S -^ln2~C / (-l) + O0r 1 ; 



Sx=i = ^ 2 + ^ln S -^ln2-3C'(-l) + 0( S - 2 ), 

S x=2 = s 2 -s + ~lns+^ln2- ^'(-1) + O^T 1 ) (A.li; 



where 



7rp R. (A. 12) 



Here we used A = (3/2 = 1/2, 1,2 instead of (3. Redefining s as in f]A.7|) and 
using (1A.8j) . we can summarize the first three terms of (1A.11I) in a compact 



form in flATTUl) . 



A. 4 Free fermions on the lattice 

For non-interacting one-dimensional fermions on the lattice (and the corre- 
sponding isotropic XY spin chain) the asymptotic behavior of EFP was de- 
rived in [24J using Widom's theorem on the asymptotic behavior of Toeplitz 
determinants. Introducing the Fermi momentum kp = irpo and using units in 
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which the lattice spacing is 1 we have 

(^ln2 + 3C'(-l))+0(iT 2 ). 

(A.13) 

In the continuous limit kp — > the (IA.13I) goes to its continuous version flA.2j) . 

A. 5 Bosons with delta repulsion 

The model of bosons with short range repulsion is described by 

1 N 

H =2^2 p2 i +9 S S ( x j~ x k), (A.14) 

3=1 l<j<k<N 

where g is the coupling constant. This model is integrable by Bethe Ansatz 
|63j . It was derived (conjectured) in [M] that the leading term of the EFP is 

S=±(KR)*[1 + I(g/K)}, (A.15) 

where K is the Fermi momentum in the Lieb-Liniger solution [S3] and 

I(x) = ^ L 7^7 L vt^ log W +<»-*)»)■ (A - 16) 

The limit /(x — * oo) = corresponds to the free fermion result ( 1A.2I) 
(Tonks-Girardeau gas), while the limit /(x — > 0) = 1 is the result for the EFP 
of free bosons. 



S = -4R 2 In cos — + - In 
2 4 



2i?sin^ 
2 
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A.6 XXZ model 

The Hamiltonian of the XXZ model is given by 

h = -jH >';-:;. i + i + , (a. I?) 

where the sum is taken over the sites j of a one-dimensional lattice and a x,y,z 
are Pauli matrices. 

Let us parametrize the anisotropy as 

A = COS7HA (A. 18) 

The asymptotic behavior of the EFP as n = 2R — ► oo was found in [281 129] 

P{n) ~ A^C"™ 2 , (A.19) 



ix\f^2/K I Jo * cosh(2tz/) sinh(t) r 



where 



and the exponent 7 was conjectured in [29] to be 



v 1 



7 = — + — -. (A.21) 

' 12 3(l-i/) V ' 
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Appendix B 

Asymptotic behavior of Toeplitz 
Determinants 

As we showed in Chapter [3[ the asymptotic behavior of the EFP for ( 12. ip at 
n — > oo is exactly related to the asymptotic behavior of the determinant of the 
corresponding Toeplitz matrix ( 13.8113.9113. 101) and can be extracted from known 
theorems and conjectures in the theory of Toeplitz matrices. These types of 
calculations have been done first in [321 [31] for spin-spin correlation functions. 
It is well known that the asymptotic behavior of the determinant of a Toeplitz 
matrix as the size of the matrix goes to infinity strongly depends upon the 
zeros and singularities of the generating function of the matrix. 

A very good report on the subject has been recently compiled by T. 
Ehrhardt [65]. Here we want to recapitulate what is known about the de- 
terminant 

D n [a] = det(S n ) = det \s(j - k)\ n hk=Q (B.l) 
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ofan + lxn+1 Toeplitz matrix 



/s(0) s(-l) s(-2) 
*(1) s (0) 
*(2) *(1) *(0) 



s(-n) \ 
s(l-iV) 
s(2 - iV) 



(B.2) 



\s(ra) s(n-l) s(n-2) ... s(0) / 



with entries generated by a function o~(q): 



S (l) 



.(,)e-g, 



(B.3) 



where the generating function <r(g) is a periodic (complex) function, i.e. a(q) = 
a{2ix + q). 

In this work we dealt only with generating functions with zero winding 
number 

Inda(g)= f p- A log a(q) = (B.4) 

and this brief review will be limited to this condition. This was not the case 
in the study of Barouch et al. [51] . where the generating function (see foot- 
note after ( 13. 1011 ) had non-zero winding number in some regions of the phase 
diagram. 



B.l The Strong Szego Theorem 

If a(q) is sufficiently smooth, non-zero and satisfies Inda(g) = (i.e., the 
winding number is 0), we can apply what is known as the Strong Szego Limit 
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Theorem ([66], [6?]), which states that the determinant has a simple exponen- 
tial asymptotic form 



D n [a] ~ E[a]G[a} n n -> oo, (B.5) 
where G[cr] and E[o~] are defined by 

oo 

G[a] = exp (5"o, -E^cr] = exp fccr fc o-_ fc (B-6) 

fe=i 

and <jfc are the Fourier coefficients of the expansion of the logarithm of cr(g): 

oo 

\oga{q) = °kd kq . (B.7) 



k=— oo 



B.2 The Fisher-Hartwig Conjecture 

Over the years, the Szego Theorem has been extended to consider broader 
classes of generating functions by relaxing the continuity conditions which 
define a "smooth function", but it remained limited to never- vanishing func- 
tions. Therefore, some extensions have been proposed to the Szego Theorem 
in order to relax this latter hypothesis. When the generating function has 
only pointwise singularities (or zeros), there exists a conjecture known as the 
Fisher-Hartwig Conjecture (FH) [68]. \_ 

When o~(q) has R singularities at q — 9 r (r = 1..R), we decompose it as 



1 This conjecture is still not completely proven. For details and status of the conjecture 
seeRef. [37|. 
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follows: 

R 

a(q) = T {q)Y[e iKr[iq - er)mod2n - n] (2 - 2cos(g - 6 r )) Xr (B.8) 

r=l 

so that r{q) is a smooth function satisfying the conditions stated in the previ- 
ous section. Then according to FH the asymptotic formula for the determinant 
takes the form 

D n [a] ~ E [r, {«„}, {AJ, {6 a }] n^-^) G [r] n n -> oo, (B.9) 
where the constant pref actor is conjectured to be 

E [t, K}, {AJ, {0J] = £[r] n^( eie ")"^ Ar ^( e " ier ) Kr_Ar 

r=l 

x ^ _ e i(6» 3 -e r )^ (Kr+Ar)(Ks " As) 

l<r^s<i? 

T"T G(l + K r + A r )G(l - K r + A r ) 

x n — Giisj — <R10) 

£/[r] and G[t] are defined as in (IB. 61) and r± are defined by decomposition 

r(g) = t_ (e") G[r]r+ (e"") , (B.ll) 

so that r + (r_) are analytic and non-zero inside (outside) the unit circle on 
which r is defined and satisfy the boundary conditions r + (0) = t_(oo) = 1. G 
is the Barnes G-function, an analytic entire function defined as 

G(z + 1) = ( 27 r) z / 2 e-[ z+ ^ +1 ^ 2 ]/ 2 II ( X + -) ~ e " 2+f ^ (B.12) 
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where 7^ ~ 0.57721 ... is the Euler-Mascheroni Constant. 

This conjecture is actually proven for some ranges of parameters n r and A r 

or fully for the case of a single singularity (R— 1), see [691 Ell- 
in many simple cases it is possible to find the factorization of r into the 

product of r + and r_ by inspection. More complicated examples like the ones 

presented in this work require a special technique to obtain this factorization, 

which is known as the Wiener-Hopf decomposition: 



logr+O) = j 



dz logr(,z) 
27ri z — w 



\w\ < 1 



logr_H = -^^) 
27ri z — w 



where the integral is taken over the unit circle. 

In light of these formulas, it is useful to present the parametrization (IB. 81) 
in a form which makes the analytical structure more apparent. Changing the 
variable dependence from q to z = e iq , we write 



a(z) = t[z, 



n >-7 fi-?) • p.") 



r=l x 



where z r = e . 



B.3 The Generalized Fisher-Hartwig Conjec- 
ture 

Despite the considerable success of the Fisher-Hartwig Conjecture, few exam- 
ples have been reported in the mathematical literature that do not fit this 
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result. These examples share the characteristics that inequivalent represen- 
tations of the form ( IB. 81) exist for the generating function a(q). Although 



no theorem has been proven concerning these generalization of the 

Fisher-Hartwig Conjecture (gFH) has been suggested by Basor and Tracy [37] 
that has no counter-examples yet. 

If more than one parametrization of the kind (IB. 81) exists, we write them 
all as 

R 

a(q) = r\q)]Je iK ^- 9r)mod2n - n] (2 - 2cos(q - 6 r )) K , (B.15) 

r=l 

where the index % labels different parametrizations (for R > 1 there can be 
only a countable number of different parametrizations of this kind). Then the 
asymptotic formula for the determinant is 

D>]~^£[^{<},{A;},{0J] n^G[rT n^oo, (B.16) 



ier 



where 



n(<) = E((a;) 2 -«) 2 ), (b.17) 

r=l 

T = ji || Re[fi(z)] = maxRe[fi(j)]| . (B.18) 



The generalization essentially gives the asymptotics of the Toeplitz de- 
terminant as a sum of (FH) asymptotics calculated separately for different 
leading (see Eq. (IB. 181) ) representations (IB. 150 . In Sec. 13.3.2.11 we used the 
sum of all (not necessarily leading) representations and showed that it also 
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correctly produces the first subleading corrections to the asymptotics of our 
Toeplitz determinant. 



B.4 Widom's Theorem 

If a(q) is supported only in the interval a < q < 2n — a as in our model for 
7 = 0, singularities are no longer pointwise and one should apply Widom's 
Theorem [38]. It states that the asymptotic behavior of the determinant in 
this case is 

D n [a] ~ 2 1 ^ e 3C'(-D (sin|)" V4 E\p]WG\p] n (co.^f , (B.19) 
where E and G are defined in (IB.6P and 

p(q) = a (2* cos -1 cos ^ cos q j (B.20) 
with the convention < cos -1 x < n. 



For the case considered in Section 13751 the generating function is constant, 
E[p] — G[p] — 1, and (IB. 191) simplifies considerably giving 

D n [a] ~ 2 1 /iV^'(- 1 ) (sinl)" 17 ^- 1 / 4 (cos|) n2 . (B.21) 
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Appendix C 

A brief introduction to the 
Bethe Ansatz 

The Bethe Ansatz technique is a very powerful tool to extract information 
about an integrable system. Its main advantage is to convert a quantum 
many body problem (the problem of solving a system of coupled Schrodinger 
equations) into the solution of a much simpler algebraic system. This solution 
specifies the energy eigenf unctions of the system, although implicitly. Even 
more importantly, it allows for the calculation of thermodynamic quantities in 
the thermodynamic limit. 

Integrable models are theories with a high degree of symmetry which guar- 
antee as many conserved quantities as the number of degrees of freedom of 
the system. Most of the results obtained in theoretical physics are derived as 
the consequence of some approximation, since the whole theory is normally 
too complicated to be solved directly. In one dimension all the systems are 
strongly interacting because of the limited dimensionality and for this reason 
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is it hard to control perturbative calculations. 

This is what makes integrable models so important in one dimension. The 
Bethe Ansatz is one of the main tools in the theory of quantum integrable 
models. It has a rich mathematical structure and a wide range of applications. 
Clearly, we cannot recount its full scope and cite the many important works 
in this subject in the space we have. We refer the reader to the recent book 
by Sutherland [7T] for its clarity and for its ability to introduce even the 
inexperienced reader to quite advanced topics. 

Here we summarize basic formulas of the Bethe Ansatz which are needed 
in the hydrodynamic approach. 

C.l The Bethe Wavefunction 

An integrable model is an exactly solvable theory. Beside this intuitive defini- 
tion, it is quite complicated to identify and formalize the exact meaning of the 
notion of integrability. We are not going to pursue the problem of identifying 
the integrability of a model and we will assume that somehow we know we are 
dealing with such a theory. 

Let us consider a one-dimensional system of N identical particles interact- 
ing with a pair potential. The Hamiltonian for the quantum problem can be 
written as: 



We assume that the potential is such that the model is integrable and we will 
not specify it further at the moment. 
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The first ingredient in the Bethe Ansatz is to assume a plane waves super- 
position form for the wavefunctions of the system: 

= A(P)e^i p h*< (C.2) 
p 

where the sum is carried over all the permutations P of the momenta of each 
particle in the system. This formula is valid when the particles are in order, 
i.e. Xi < Xi+i for every i, and allows for an easy extension to the other 
configurations by employing the symmetry or antisymmetry properties of the 
wavef unction. 

One of the fundamental conditions for a model to be integrable is that there 
are no two-particle irreducible interactions, i.e. all processes in the system 
can be viewed as a sequence of two particle scattering, and the order of the 
scatterings does not matter. This latter condition is formally expressed as the 
Yang-Baxter equation. 

The important quantity in a two-particle scattering is the phase acquired 
by the particles after the interaction. In the Bethe Ansatz construction the 
"scattering phase" 6±(k) plays a very important role. We can define a different 
phase for bosons 0+(k) and for fermions 9-{k) to take in account the different 
statistics under a particle exchange. One can prove that 9±{k) is a real odd 
function of k {0*±(k) = 9±(k*), 9±(-k) = -9±(k)), where k = k x - k 2 is the 
difference between the momenta of the two interacting particles. 

If the particles were distinguishable, we would describe the scattering pro- 
cess through the transmission and reflection amplitudes. These can be ex- 
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pressed in terms of the scattering phase as 

T(k) . , 

fl(fc) = ± . (C.3) 

It is straightforward to solve the one- dimensional Schrodinger problem for a 
given potential and derive from it the scattering phase. We skip these deriva- 
tions here (see, for instance, [71]) and just list the results for a couple of 
integrable potential we will use in our hydrodynamic approach. 

First, let us consider a contact interaction: 

U(r) = c8(r). (C.4) 

For fermions, the antisymmetry of the wavefunction prevents the particles 
from interacting, so 

0_(Jfc) = O. ( c - 5 ) 
For Bosons the solution of the scattering problem brings 

8+(k) = -2arctan(A;/c). (C.6) 

If the particles are distinguishable, one can then compute the transmission 
and reflection amplitudes from (IC.3I) as 
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Let us consider "Calogero- Sutherland" particles, i.e. particles interacting 
with the potential: 

u(r) = (c . 8) 

where A quantifies the interaction strength and the statistics of the particles. 
In fact, the scattering phase is 

6{k) = 7r(A — 1) sgn (A;). (C.9) 

One sees that A = 1 corresponds to a free Hamiltonian with antisymmetric 
wavefunction and therefore describes free fermions, while for A = the po- 
tential vanishes again but the wavefunction is symmetric and it describes free 



bosons. In analogy to higher dimensional cases, Ca 



ogero-Sutherland particles 







can be viewed as anyons with statistics given by 9, 

Once one knows the scattering phase 0(k), it is a just a tedious applica- 
tion of algebra and combinatorics to write the explicit form of the amplitudes 
A(P) in ( 1C.2I) . We are not going to pursue the problem of writing down the 
wavefunction explicitly here, especially since most of the important physical 
results can be obtained already without this knowledge. For a very detailed 
derivation and explanation of these details, we refer the reader to 



^dn fact, in one dimension there is no way to bring a particle pass another avoiding the 
interaction, so the concept of fractional statistics, introduced in two-dimensional physics, is 
of a different nature in one dimension. 
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C.2 Periodic boundary conditions 

We established that every scattering can be decomposed in a sequence of two- 
particle scatterings and that the effect of such a process on the wavefunction is 
just an additional phase that depends on the difference between the momenta 
of the two interacting particles. 

We are now in position to establish the fundamental equations and results 
for the Bethe Ansatz. Let us consider a system of N particles in a box of size 
L. Eventually, we want to take the thermodynamic limit for the size of the 
box L and the number of particles N to go to infinity in such a way that the 
density N/L stays finite. So it doesn't matter the kind of boundary conditions 
we impose and we choose periodic boundary condition on the wavefunction 
(Q: 



x 2 ,...,Xj + L,...,x N ) = ^(xi, x 2 ,...,Xj,..., x N ). (C.10) 

For a particle j to wind around the box, it would have to scatter through 
every other particle of the system acquiring this way a phase 

N 

Y[e~ ie ^- ki \ (C.ll) 

i=i 

where we remembered that 8(0) = 0. Moreover, during the motion the particle 
will acquire a dynamical phase 

e iLk i /n . (C.12) 
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Putting both effects together, periodic boundary conditions amount to impose 

N 

1 = e iLk > /h Y[e- i9{k >- kl) (C.13) 

i=i 

for each of the N kj momenta in the system. 

Taking the logarithm of ( 1C.13I) we arrive at a system of N coupled algebraic 
equations 

N 

2irlj = Lkj/h + e ( k j ~ h), (C.14) 

i=l 

which are called the "Bethe Equations" or the "Fundamental Equations". The 
N integers Ij are the winding numbers of the phase and are effectively the 
quantum numbers for the state described by the wavefunction. 

The importance of the equations (1C.14|) is paramount, since they translate 
the problem of solving a complicated differential equation into a system of 
algebraic equations. Their solution gives the momenta kj for the system and 
this information is already sufficient to calculate the total momentum 

N 

P = J> (C.15) 

and energy 




(C.16) 



i=i 



It is important to note now that, although our derivation would lead to 
the conclusions that the kj can be identified as the actual momenta of the 
individual particles, this interpretation is in general not correct. Not only is 
it incorrect to assign a defined momentum to a particle in a region where the 
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particle is interacting, but, more fundamentally, more complicated models will 
assign quasi-momenta labels kj to their degrees of freedom, but these labels 
cannot be interpreted as momenta in a traditional way. 

In general, the quasi-momenta kj should be viewed as points in an appro- 
priate phase-space that describe the system. The integers Ij will specify the 
state of the system and it can be shown that these N integers have to be 
chosen all distinct (if this was not the case, two particles would have the same 
momentum and would not interact and this would make the system singular). 

C.3 Zero temperature thermodynamics 

The set of integers Ij in (IC.14I) defines the state described by the wavefunction. 
The ground state, the state with the lowest energy, can be identified with the 
state with the set of Ij running from —N/2 to N/2, N being the number of 
particles in the system. Clearly, from flC.15j) we see that this state has zero 
momentum. By choosing a different set of integers, one constructs the excited 
states of the model. 

When we take the thermodynamic limit, we let the number of particles N 
and the length of the system L go to infinity as we keep their ratio, the density 
of particles, finite 

N/L = p. (C.17) 

As we increase N, one can show that the distribution of the momenta kj of 
a system grows more dense. In the thermodynamic limit one can prove that 
the average distance between two neighboring momenta scales like 1/L and 
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we can define the distribution function 



r(k j )= lim — — >0 (C.18) 

N— >oo,L^oo L{Kj+i — Kjj 



which defines the density of particle in this quasi-momentum space. 
We can write the Bethe Equations ( 1C.14I) as 

1 - 

i=i 

where we defined the "counting function" y(kj) = which is a monoton- 
ically increasing function that counts the integers as a function of the quasi- 
momenta. By definition, 



yikj) - y(h) = ^(Ij - Ii) (C.20) 



and one can show that 



*r(*) = ^, (0.21) 



establishing a direct connection between the distribution of the integers and 
of the quasi-momenta. 

In the thermodynamic limit, the system of algebraic equations ( 1C.14j) can 



be written as an integral equation for the counting function and the momentum 
distribution: 

y(k) = k+ [ ^ 6{k - k')T(k')dk' (C.22) 



174 



and, by taking the derivative of this equation by k, 



1 1 pttmax 

r(k) = _ + _/ 6'{k - k')r(k')dk' 

Z7l 2ll h . 

,J "•rain 
1 pkmax 

= — + / K{k - k')T(k')dk' (C.23) 
2tt II, . 

where we introduced the kernel of the integral equation as the derivative of 
the scattering phase: 

Equation (1C.23I) allows us to determine the distribution of the quasi-momenta. 
This distribution depends on the support of the kernel, in equation ( 1C.23I) the 



limits of integration fc min and k max . The support is determined by the choice of 
the integers in the original equations (10. 14[) . For the ground state, the limits 
of integration are symmetric ik m in = —k ma x = kp)- A direct way to determine 
the limits of integration is to calculate the number of particles per unit length: 

N/L = p= f F r(k)dk (C.25) 

and invert this equation to calculate kf in terms of N. 

Finally, we can write AO. 15D in the thermodynamic limit as 



/kp 
k r(k)dk = 0, (C.26) 
■k F 



where we have used the fact that r(—k) = r(k), and we rewrite (1C.16I) as 



"kp .2 

E/L= I —r(k)dk. (C.27) 

-k F 2 
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Equations (1C.25I) and (1C.27I) . together with (1C23j) . define the equation of 
state of the system: 

pe(p) = j. (C.28) 
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Appendix D 

Integrability of Gradient-less 
Hydrodynamics 

In this appendix we construct the conserved quantities of gradient-less hy- 
drodynamic theories. It is a little-known fact that if we consider a Galilean 
invariant system and we neglect terms containing spatial derivatives of the 
density or of the velocity, we can construct an infinite series of conserved 
quantities ( "integrals of motion"). 

Although this result has been derived in the 1980's in the study of the 
mathematical structure of a class of differential equations (defined as of "hy- 
drodynamic type" [76]), it is not well-known in the physics community Math- 
ematicians have traced the source of these conserved quantities to a property 
known as "multi- Hamiltonian structure" [73]- [75]. This means that there ex- 
ists more than one set of Hamiltonian and symplectic structure, i.e. Poisson 
brackets, that generates the dynamical equations of the model. 

Besides the Hamiltonian (14.51) which in connection with the Poisson brackets 
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(14. 8 p generates the equations of motion (14.9114. 101) . we can construct a differ- 
ent Hamiltonian and Poisson brackets that would lead to the same dynamical 
equations. This interesting property allows integrals of motion to be trans- 
lated between the different Hamiltonian systems and this generates a ladder 
structure in which one can generate new conserved quantities using the ones 
of the other system and vice versa. 

We are not going to describe this structure any further, since what is 
missing in our opinion is a clear physical understanding of the origin of all 
these conserved quantities. We know that in general conserved charges are 
due to some symmetry of the model, but we have not been able to identify 
what symmetry would guarantee such an infinite series, actually, a double 
infinite serieJl]. 

We know that the free fermions model possesses a (or even Gl(oo)) 
symmetry and this is the source of the integrability of the system. We think 
that even interacting hydrodynamic theories could possess the same vast sym- 
metry and that, in some sense, by neglecting gradient corrections we dis- 
card most of the content of the system and end up with some essentially free 
fermions in disguise. It might be possible to establish an exact mapping be- 
tween any two hydrodynamic theories (and in particular to map an interacting 
model into the free fermion model), but so far we failed in constructing such 
transformation. 

After introducing the model in section \D.1\ we discuss the integrability of 

the system in section ID. 21 and argue that the presence of an infinite number 

1 Certainly, any spatial symmetry like Galilean invariance is not sufficient to guaranty 
such a high number of conserved quantities, since it is not infinite-dimensional. 
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of conserved densities does not guarantee the full integrability of the theory. 
In section ID. 31 we explicitly construct the double infinite series of integrals 
of motion. As far as we know, this is the first time that such an explicit 
construction has been reported. Finally, we are going to discuss these results 
in section lD.41 

D.l Hydrodynamic Hamiltonian and equations 
of motion 

Let us consider a generic Hamiltonian of the Hydrodynamic type like the ones 
consider in Chapter |H We require the theory to have Galilean invariance and 
neglect terms with derivatives of the fluid density p or of the fluid velocity v: 

n = ^ + P e(p), (D.l) 

where e(p) is the internal energy of the fluid. The theory is defined by spec- 
ifying the Poisson brackets between the density p and the velocity v as in 
AMD: 

{p{x),v(y)} = -i6'{x-y), (D.2) 

while the other Poisson brackets vanish identically. 

From the Hamiltonian and the Poisson brackets we find the equations of 
motion to be 

Pt = -pv x - vpx = -{pv) x , (D-3) 
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v t = -vv x - (pe) pp p x = -vv x -p x , 



(DA) 



where 



v 




pp 



(D.5) 



is the "Sound Velocity" of the fluid and where we adopted a notation for which 
8aB = Ba- Eq. (1D.3P expresses the continuity condition for the fluid, while 
Eq. f]D.4|) is the Euler dynamical equation of motion. 

D.2 Integrability of the Hydrodynamic theory 

Some integrals of motion of the theory are trivial: p, pv and H are obviously 
conserved densities. 

To find more conserved quantities let us start with 



and find under which conditions is the function f(p, v) a conserved density. To 
this end, we calculate the commutation relation between / and H = J Tidx 
and we impose it to be zero: 




(D.6) 




(D.7) 
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where g{p,v) is some unknown function. Equating the second and fourth 
line and imposing the condition g pv = g vp , we can eliminate g{p,v) from the 
equations and find that f(p, v) is a conserved density if 

v 2 s 

fp P = -4/w (D-8) 
P 2 

Eq. (ID. 81) allows us to find integrals of motion. To determine how many 
of them can be simultaneously specified, let us calculate the commutation 
between two conserved integrals 

h = Jdxf(p,v) I 2 = Jdxh(p,v) (D.9) 

with conserved densities f(p,v) and h(p,v) satisfying Eq. ( ID. 81) : 

{h,I 2 } = [ dxd x g(p,v)=0. (D.IO) 



Solving this condition as before we find that the two integrals of motion com- 
mute if 



which is identically satisfied since both f(p,v) and h(p,v) satisfy Eq. (1D.8j) . 



Therefore, we see that any solution of ( 1D.8|) is an integral of motion and, 
since we can construct an infinite series of solution, linearly independent from 
each others, the theory admits an infinite series of mutually commuting inte- 
grals of motion and their respective conserved densities. One would conclude 
from this that the theory is integrable, but it is not necessarily the case. 
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To be more precise, there are many definitions of integrability for a system 
with an infinite number of degrees of freedom. If by integrable we mean that 
every degree of freedom admits a representation in terms of action-angle vari- 
ables (integrability according to Liouville), the existence of an infinite number 
of integrals of motion (the "actions" conjugated to the "angle" variables) might 
not be enough to exhaust all the degrees of freedom. 

Therefore, whether gradient-less hydrodynamic theories are integrable ac- 
cording to Liouville or only integrable in the lesser meaning of possessing 
infinitely many conserved quantities is still not clear. 

D.3 The integrals of motion 

To understand better the structure of the integrals of motion, let us start with 
the simplest ansatz for a solution of Eq. ( ID. 81) : 




(D.12) 



so that Eq. ( ID. 81) becomes 




(D.13) 



Expanding this solution in a Taylor series in powers of k: 




(D.14) 



n 



we can calculate 




(D.15) 
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Using Eq. (ID. 121) we can calculate the same quantity as 



9 v f(p,v) 



Ke~ KV J K (p) 



(D.16) 



n 



and we find that 



f n (p,v) = d v f n+1 (p,v). 



(D.17) 



Since f(p,v) is a conserved density, so are the coefficients of the Taylor 
expansion f n (p,v). We showed before that any integral of motion commutes 
with any other; therefore, the functions f n (p,v) are the conserved densities of 
a series of mutually commuting integrals of motion. 

The choice of the ansatz in Eq. (ID. 12h provides only two linearly inde- 
pendent integrals of motion, but by expanding any of these two solutions in 
a Taylor Series, we recover an infinite series of conserved densities. It can be 
shown that there exist exactly two infinite series of integrals of motion and 
each of them can be generated from the two solutions of Eq. (ID. 121) . Moreover, 
we found a recurrence relation among each series, Eq. ( ID. 171) . We can use this 
relation to generate the entire series from the first element in the series (the 
one with the lowest power in v) and successively integrating in v; we only need 
to determine the integration constant. 

Let us calculate the elements of the first series. We already know three of 
them: 



-i 
o 




(D.18) 
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dxpv 



dx 



(D.19) 
(D.20) 



The general structure of this series is 



n j ^ X 3n 



(P. V ) 



(D.21) 



with 



■ i pv 

Jn 

n 



[n/2] 



n-2k 



k=l 



(D.22) 



where [n/2] means the highest integer smaller or equal to n/2. 

We can determine the coefficients 0^(p) by requiring dtj\ = d x g for some 
function g(p,v). Solving this condition brings: 



0n = (n-l)pe(p) 



(D.23) 



- > n)pp 



(n - 2A; + 2)(n - 2A; + l)^^" 1 



(n-2A; + 2)(n-2A; + 1) 



n 2 ^™ 



fc = 2 . . . [n/2] 



(D.24) 



and in this way we can recurrently generate the series. For instance: 



dx 



pv 



2pe(p)v 



(D.25) 



It is easy to show that the velocity v is also a conserved density and it 
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constitutes the first element of the second series of integrals of motion: 



(D.26) 



with 



[n/2] 



Jn 



n 



,n-2k 



k=l 



We can determine the coefficient as before and find 



fn- 1 



/>/> 



P 



(n-1) 



P z 



{n-2k + 2)(n - 2A; + l)^^" 1 



and generate the whole series: 



(n-2/c + 2)(n-2/e + l)%^- 1 

P 



(D.27) 



(D.28) 



fc = 2... [n/2] 



(D.29) 



i 3 



dxu 
die 
dx 



r v 2 



£(p) 



(D.30) 
(D.31) 
(D.32) 



where 



(D.33) 

It is interesting to notice that the conserved densities j l n , with the given 
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definition for the coefficients, naturally follow the condition in Eq. (1D.17|) . i.e. 



j l n (p,v) = d v j l n+1 (p,v). (D.34) 

and it can be shown that they are the coefficients of the Taylor expansion of 
the two solutions of Eq. (ID.12I) . 

D.4 Conclusions and open questions 

We have shown that a hydrodynamic theory having Galilean invariance and 
without gradient corrections has two infinite series of integrals of motion. We 
still do not understand completely the nature of this integrability, i.e. the 
symmetry group it comes from. 

The hydrodynamic description of free fermions has the Hamiltonian: 



H = + ^ (D.35) 
2 6 H K J 



with sound velocity 



v s = up. (D.36) 



It is known to be integrable, its symmetry group being W^, and its integrals 
of motion being simply: 

J dx^(v±p) k (D.37) 

which are linear combinations of the l\ and 1% (1D.2HD.26I) . 

It is natural to think that the origin of the integrability of the generic 
Hamiltonian in Eq. (ID.lj) is the same as for the free fermions, but this has not 
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been confirmed yet. 

It is worth noticing, however, that unlike other integrable theories this 
system does not have solitons and that any wave eventually develops singular- 
itiej^l. This is in sharp contrast with the intuitive concept of integrability and 
is probably a sign that the double infinite series of integrals of motion is not 
sufficient to completely describe and constrain the dynamics of the system. 

A lot of work has been done on equations of the hydrodynamic type and a 
good understanding of the integrability of these theories in terms of an under- 
lying multi-Hamiltonian structure is established (ref. [73]- [75] for a very short 
and non-comprehensive list), but we believe that the physical nature of the 
integrability, i.e. the underlying symmetry group, is not yet been explained. 
Reference [76] has the most symmetry-oriented approach. 

It is also known that there exist integrable hydrodynamic theories with gra- 
dient corrections, namely the ones describing the Calogero-Sutherland model 
and the free bosons with delta repulsion (Lieb-Liniger bosons). It would be in- 
teresting to understand what preserves the integrability of this theories, which 
were integrable even before the addition of the gradient corrections, as we just 
showed. 



2 Systems like these can be viewed as "dispersionless limits" of other integrable theories 
and can therefore be considered singular, since more than one integrable model can have 
the same dispersionless limit. 
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